Dynamics of Package Cushioning 

By RAYMOND D. MINDLIN 

Introduction 

MECHANICAL damage is a common occurrence in the transportation 
of packaged articles. The causes of failures are generally inadequate 
protective cushioning, lack of ruggedness of the outer packing container, 
or occasional abnormal weakness of the packaged article. The first of these 
difficulties is the subject of this paper. 

One of the major influences in reducing the incidence of mechanical failures 
of packaged articles in recent years has been the use of the drop test. The 
drop test is performed simply by raising the package to a specified height and 
dropping it to the floor. The package and its contents are then examined 
for damage. This is a go-no-go test and requires a large number of samples 
before a reliable estimate of quality can be made. An adequate number of 
tests is prohibitive when the article packaged is costly. In such cases it is 
important, and in any case it is useful, to supplement the drop test data with 
measurements and calculations. It is also possible to evolve rational pro- 
cedures for designing packages, as described in the present paper, so that 
a particular product will survive a drop test at any specified height, with a 
known factor of safety and with a minimum amount of space assigned for 
cushioning. The drop test then becomes only a check instead of playing an 
integral role in a cut and try design procedure. 

Assuming that the outer container is adequate, the survival of a packaged 
article in a drop test still depends upon a large number of factors descriptive 
of the mechanical properties of both the cushioning medium and the pack- 
aged item. However, the more important properties can be grouped so 
that they may be replaced by knowledge of only the following factors: 

(1) The magnitude of the maximum acceleration that the cushioning 
permits the packaged item to reach. 

(2) The form of the acceleration-time relation. 

(3) The strengths, natural frequencies of vibration and damping of the 
structural elements of the packaged article. 

Part I of this paper is concerned primarily with methods for predicting 
maximum acceleration of the packaged article with emphasis on non-linear 
cushioning. Part II deals primarily with the prediction of the form of the 
acceleration-time relation. Part III deals with the effect of acceleration on 
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the packaged article and gives methods for determining whether or not the 
strength of the packaged article will be exceeded. The strength deter- 
minations themselves are not dealt with here; but the information in Part III 
is essential in interpreting and applying the data obtained in strength meas- 
urements. In Part IV some consideration is given to the influence of dis- 
tributed mass and elasticity. 

It cannot be emphasized too strongly that the determination of the mech- 
anical properties of the packaged article, not dealt with in this paper, is an 
essential preliminary to a rational design procedure for packaging. The 
whole purpose in designing package cushioning is to limit the forces which 
may act on the packaged item. If one does not know to what values to 
limit the forces, a rational design procedure cannot be applied. 

It is interesting to observe that the methods described here for analyzing 
and designing package cushioning are directly applicable to the design of 
shock mounts intended to protect equipment from the effects of a sudden 
change in velocity. All of the principles, formulas and design curves given 
here may be used in the shock mount problem with the simple substitution 
of V 2 /2g for h, where h is the height of drop in the packaging problem, g is 
the acceleration of gravity and V is the velocity change in the shock mount 
problem. 

This paper is essentially a report on a study undertaken at the Bell Tele- 
phone Laboratories, Inc., in the Electronic Apparatus Development 
Department. The results have been applied to the packaging of large 
vacuum tubes and all of the examples used to illustrate the analysis and 
design procedures in the paper are taken from vacuum tube applications. 

Miss H. A. Lefkowitz, Member of the Technical Staff, Bell Telephone 
Laboratories, assisted in the mathematical studies. The oscillograms, used 
as illustrations, were prepared under the supervision of Mr. F. W. Stubner, 
Member of the Technical Staff, Bell Telephone Laboratories. Figure 
3.8.2 was taken from a thesis submitted by Mr. C. Ulucay in partial fulfill- 
ment of the requirements for the degree of Master of Science in the Depart- 
ment of Civil Engineering at Columbia University. The calculations for 
Figs. 3.5.1 to 3.5.6 and Fig. 3.2.2 for fii > were performed on the Westing- 
house Mechanical Transients Analyzer under the supervision of Dr. G. 
D. McCann, Transmission Engineer, Westinghouse Electric and Manu- 
facturing Company. 

Assumptions 

The procedures to be described for the analysis and design of package 
cushioning are based on applications of a few simple laws of mechanics to 
an idealized mechanical system representing the package and its contents. 



DYNAMICS OF PACKAGE CUSHIONING 



355 



Essentially, a package consists of 

1. Elements of the packaged article which are susceptible to mechanical 
damage. 

2a. The packaged article as a whole. 

2b. A cushioning medium (excelsior, cardboard spring pads, metal springs, 
etc.) 

3. An outer container (cardboard carton, wood packing case, etc.) 
The four major components are illustrated schematically in Fig. 0.2.1. 
The system is further idealized by "lumping the parameters"; for example, 
the outer container is considered as a single mass, the cushioning is con- 
sidered as a massless spring with friction losses. The result of this idealiza- 
tion is to lose some of the fine detail of the real distributed system such as 
wave propagation through the cushioning and higher modes of vibration in 
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Fig. 0.2.1 — Schematic representation of a package. 
1. Element of packaged article 
2a. Packaged article as a whole 
2b. Cushioning 
3. Outer container 

the package structure and in the packaged article. Some consideration of 
these details is given in Part IV. 

The idealized system is illustrated in Fig. 0.2.2. The major components 
of the system are as follows: 

1. A structural element of the packaged item is represented by a mass 
(mi) supported by a linear massless spring with or without velocity 
damping. The mass m\ is assumed to be small in comparison with 
the mass of the whole packaged item. 
2a. The whole packaged item is represented by a mass mi. 
2b. The cushioning is represented by a spring which may have a linear 
or non-linear load-displacement characteristic and which dissipates 
energy through velocity damping or dry friction. Permanent de- 
formation of the cushioning is not considered, that is, in a repetition 
of the drop test it is assumed that the package has the same properties 
as before the first test. A properly designed package will have essen- 
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tially this characteristic. The mass of the cushioning is assumed to 
be small in comparison with mi, except in Section 4.2. 
3. The outer container is represented by the mass m 3 . The impact of 
m 3 on the floor is assumed to be inelastic and during contact the rela- 
tive displacement between m 3 and the initial position of the floor 
is assumed to be small in comparison with the relative displacement 
between Wo and m 3 . In other words, no spring action is assigned to 
the outer container and the floor is considered rigid. 



Element of 

Packaged 

Item 

Packaged 
Item 



Cushion 



Outer 
Container 



Height 

of drop h 




///////////////////// 

(a) (b) (c) 

Fig. 0.2.2 — Idealized mechanical system representing a package in a drop test. 

PART I 
MAXIMUM ACCELERATION AND DISPLACEMENT 

1.1 Introduction 

Most of Part I is concerned with the prediction of the maximum accelera- 
tion that the cushioning permits the packaged article (m 2 ) to attain. In 
many instances this will be all the information necessary for judging the 
suitability of a cushioning system. It will be all that is necessary if the 
shape and scale of the acceleration-time function satisfy certain criteria 
which are treated in detail in Parts III and IV. If these criteria are satisfied, 
the effect of the drop on the packaged article is found by multiplying the 
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dead load stresses (obtained in the usual manner) by the ratio of the maxi- 
mum acceleration to the acceleration of gravity. If the criteria for the use 
of maximum acceleration alone are not satisfied, then Parts II and III will 
supply a numerical factor (the Amplification Factor) by which the maximum 
acceleration should be multiplied, and the remainder of the procedure is the 
same as before. 

The determination of the maximum acceleration is founded on a knowledge 
of the load-displacement characteristics of the cushioning. When the cush- 
ioning system is simple enough, the load-displacement relation may be found 
or designed by purely analytical procedures. The tension spring package, 
discussed in Sections 1.7 and 1.8, is an example where such a treatment is 
possible. In many instances, as with distributed cushioning, the load- 
displacement relation is more easily found by test. 

A load-displacement test is made by applying successively increasing 
forces, with weights or in a load testing machine, to the packaged item 
completely assembled in its package, and measuring the corresponding 
displacements. The force is applied usually by means of a rod inserted 
in a hole cut through the outer container and the cushioning to the packaged 
item. It is convenient to use a low loading rate in the test, and, in doing so, 
the effect of resisting forces that depend on velocity is lost. These forces 
are often of little importance but, in certain designs, it is necessary to con- 
sider them. This is done for velocity damping in Sections 2.5, 2.6, 3.2 and 
3.5. 

Most of Part I is concerned with cushioning having non-linear load-dis- 
placement characteristics. Linear cushioning is rarely encountered, but 
it will be treated first because of its simplicity and because it will be con- 
venient later to express the maximum acceleration in non-linear cases in 
terms of the maximum acceleration in a hypothetical linear case. 

1 .2 Derivation of Equations of Motion 

To introduce the method of analysis that will be used in Part I, the sim- 
plest possible system is considered first. The m\ system is omitted entirely, 
the mass of the outer container (w 8 ) is neglected, and the cushioning is 
assumed to have no damping or friction. There remain only the mass 
m 2 (the mass of the packaged item alone) and the supporting spring, as 
shown in Fig. 1.2.1. If the spring is linear its displacement is proportional 
to the applied load throughout the range of use (see Fig. 1.4.1). The spring 
rate (£ 2 ) of a linear spring is a constant usually expressed in terms of pounds 
per inch. The force (P) transmitted through a linear spring is therefore 
given by 

P = k 2 x«, (1.2.1) 
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where x 2 is the displacement of m 2 measured downward from its position at 
first contact of the spring with the floor (see Fig. 1.2.1). For a non-linear 
spring P will be some other function of x 2 : 



P = P(x 2 ). 



(1.2.2) 



To write the equation of motion for the mass m 2 , we consider the forces 
acting on it at any instant. These are (see Fig. 1.2.2(b)) the spring force 
P and the weight m 2 g, where g is the acceleration of gravity. When #2 
is positive (i.e., a downward displacement of m 2 from its position at first 
contact of the spring with the floor) the spring exerts an upward force P 
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(a) (b) (c) 

Fig. 1.2.1 — Elementary system. 





(a) 



(b) 



Fig. 1.2.2 — Free body diagram for elementary system. 

(a) Spring not in contact with floor. 

(b) Spring in contact with floor. 

on the mass, opposing the weight. The total downward force on m 2 is 
thus m 2 g — P. By the second law of motion, the product of the mass and 
its acceleration at any instant is equal to the applied force: 



m 2 x 2 = m 2 g — P, 



(1.2.3) 



where the symbol x 2 , representing the acceleration of m 2 , stands for the 
second derivative of displacement with respect to time (d 2 x 2 /dt 2 ). Equation 
(1.2.3) is the law governing the motion of m 2 as long as the spring is in con- 
tact with the floor. When the spring is not in contact with the floor, it can 
exert no force on the mass so that, in writing the equation of motion that 
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governs before or after contact, the free-body diagram of Fig. 1.2.2(a) should 
be used. Then 



X 2 = 



(1.2.4) 



Equation (1.2.4) holds (neglecting air resistance) from the instant the 
package starts to fall until the instant it strikes the floor and from it we 
can find the package velocity at the instant of first contact. Integrating 
(1.2.4) with respect to time, we find 

x 2 = gt + A, (1.2.5) 

where x 2 is the velocity (dx 2 /dl) and A is a constant of integration whose 
value is found from the initial condition that when / = (the instant of 
release) x 2 = 0. Thus ^4=0 and 

x, = gt. (1.2.6) 

Integrating again, 

*»-*«*+*. d-2.7) 

The value of the integration constant B is found from the initial condition 
that x 2 = — h (the height of drop) when t = 0. Hence B = —h and 

x 2 = hgt'-h. (1.2.8) 

At the instant of contact, x 2 = and, from (1.2.8), the time at first contact 
is given by tl = 2k/ g. Substituting this value of t in (1.2.5) we find, for the 
velocity at first contact, 

[*2]* 2 -o = \/2g~h- (1-2.9) 

We now have the initial conditions for finding the values of the integration 
constants in the solution of equation (1.2.3), which we proceed to obtain._ 

First multiply both sides of (1.2.3) by dx 2 /dt and write .i : 2 = - ( -^ J: 

dx 2 d (dx 2 \ . r ,dx 2 dx 2 n-nt/w 

W2-T7 t. l-rr) + -P-n- = mtK-ji (1.2.10) 



dt dt\dt ) dt dt 

or 

. d /dx 2 \ 2 . n dx 2 dx 2 

Multiplying by dt and integrating once: 

\m 2 x\ + f 2 P dx 2 = f ' m 2 g dx 2 4- C, (1.2.11) 
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where C is a constant of integration whose value is determined by the initial 
conditions that x\ = 2gh and x 2 = at the instant of contact. Hence 

f° 
C = m 2 gh -{-IP dx 2 . 

Substituting the above value of C in (1.2.11), we have 

%tn 2 xl + / P dx 2 = m 2 g(h + x 2 ). (1.2.12) 

It may be observed that (1.2.12) is an energy equation in which the 
terms have the following meanings: 

%m 2 x\ is the instantaneous kinetic energy of m 2 , 

I P dx 2 is the energy stored in the spring at any instant. It is also 
Jo 

equal to the area under the load-displacement curve up to 

the displacement x 2 , 

m 2 g(h + x 2 ) is the potential energy of the mass at its initial height h -\- x 2 

above the instantaneous position x 2 . 

Hence (1.2.12) expresses the law of conservation of energy. 

Ordinarily h is very much larger than :v 2 so that we may write, with good 

accuracy, 

\m 2 x 2 + / P dx 2 = m 2 gh. (1.2.13) 

Jo 

To the same approximation, equation (1.2.3) becomes 

m 2 x 2 + P = 0. (1.2.14) 

Equation (1.2.14) and its first integral, equation (1.2.13), are convenient 
forms for calculating events at any time during contact. Their use will be 
illustrated in Part II. For calculating only maximum displacement and 
acceleration, the equations become simpler. Let 
W 2 = weight of the packaged article ( = m 2 g), 
dm = maximum displacement of the packaged article, 
G m = absolute value of maximum acceleration of the packaged article 

in terms of number of times gravity (G m = | x 2 /g \ mayi ), 
P m = maximum force exerted on packaged article by cushioning. 
We shall limit our study to the practical regions where P > when 
a-. > 0. Then it may be seen from (1.2.13) that x 2 is a maximum when 
^2 is zero, hence 

Pdx 2 = W 2 h, (1.2.15) 



/ 

Jo 
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and, from (1.2.14), 

W 2 

where P m is the maximum value of P. If P(x 2 ) is a monotonic function, 
P m may be obtained from (1.2.2) by substituting d m for * 2 : 

P m = P(d m ). (1-2.17) 

In the unusual case where P(x 2 ) is not monotonic, the maximum value of P 
in the interval < x 2 < d m must be chosen instead of equation (1.2.17). 
The general procedure is to calculate d m from (1.2.15), P m from (1.2.17) 
and then G m from (1.2.16). If P can be expressed analytically in terms of x 2 
and if the integral in (1.2.15) can be evaluated in terms of elementary func- 
tions, simple formulas can be found for d m and G m . If this is not possible, 
then the integration can be performed graphically or numerically. Both 
of these procedures will be illustrated. In either case the maximum accel- 
eration and displacement are obtained in terms of the weight of the pack- 
aged item, the height of drop and parameters descriptive of the load-dis- 
placement characteristics of the cushioning. 

1.3 Linear Elasticity 

For cushioning with a linear load-displacement relation, equation (1.2.1) 
applies. Substituting this value of P in (1.2.15), and performing the in- 
tegration, we find 

4. = ^. (1-3.1) 

From (1.3.1) and (1.2.17), 

P m = V2hW 2 k 2 , (1-3.2) 

and, from (1.3.2) and (1.2.16), 

G„= a/W*. ' (1-3.3) 



M'2 

Notice that equation (1.3.3) holds only if there is space available for a 
displacement d„, and if the cushioning is linear and capable of transmitting 
a force P m . Also, from (1.3.3) and (1.3.1), 



and 



*. = P (1-3.4) 

Cr m 



, = M» = 2 ^. (1.3.5) 

2 2h\ dl 
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Example: Find the properties of the linear cushioning required so that 
the maximum acceleration will be 50g in a 3 ft. drop of a 20 lb. article. 
From (1.3.4), 

necessary travel, d m = — — — =1.44 inches. 

From (1.3.5), 

. , 20 X (50) 2 , n . _ , 
spring rate, k 2 = K ' = 694 lbs/in. 

From (1.2.16) 
Maximum force P m = 20 X 50 = 1000 lbs. 

1.4 Cushioning with Non-Linear Elasticity 

In practice it is rarely that a packaging system has linear spring charac- 
teristics. Departure from linearity may be due to 




■■2 

Fig. 1.4.1 — Linear elasticity. Class A. 

1. Non-linear geometry, such as in the tension spring package described 
in Section 1.7. 

2. Non-linear characteristics of distributed cushioning materials such as 
excelsior and rubber. 

3. Abrupt change of stiffness such as occurs if the packaged item can 
strike the wall of the container. 

For the purpose of developing design formulas it is desirable to have 
analytical functions to represent load-displacement characteristics. It is 
not feasible to have only one family of functions with adjustable parameters 
to fit all possible shapes of load-displacement curves. Therefore, all the 
practical shapes have been divided into six general classes, most of which 
are associated with simple functions having one or two adjustable param- 
eters. The six classes are as follows: 

Class A — Linear Elasticity. This has already been treated. Its load- 
displacement function is 

P = ta. (1.4.1) 
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Class B — Cubic Elasticity. This includes cushioning which does not bot- 
tom in the anticipated range of use, but the slope of the load-displacement 
function generally increases with increasing displacement as in the curved 
full line of Fig. 1.4.2. A suitable load-displacement function is 



p = k x 2 + rxl. 



(1.4.2) 



ko is the initial spring rate of the cushioning, as shown by the slope of the 
dashed straight line in Fig. 1.4.2, and r determines the rate of increase of the 
spring rate. The same function can be used if the slope of the curve de- 
creases gradually with increasing load as shown by the curved dashed line 
in Fig. 1.4.2. In this case the parameter r is negative. 





Fig. 1.4.2 Fi g- 14 - 3 

Fig. 1.4.2— Cubic elasticity. Class B. 
Fig. 1.4.3— Tangent elasticity. Class C. 

Class C— Tangent Elasticity. Cushioning that bottoms, but not very 
abruptly, can be represented by the load-displacement function 



_ 2k d b . irx 2 
P = tan ;r-r 

IT Idh 



(1.4.3) 



Referring to Fig. 1.4.3, h is the initial spring rate and d b is the maximum 
available displacement. The figure shows how the stiffness of the cushion- 
ing (i.e., the slope of the curve) increases as the displacement approaches 
the maximum available (d b ) at hard bottoming. The shape of the curve 
is typical of load-displacement curves for a great variety of packages with 
distributed cushioning. 

Figure 1.4.7 illustrates the wide variety of shapes of non-linear cushioning 
characteristics that can be obtained with the single function given by equa- 
tion (1.4.3) simply by varying the parameter h; and a similar set is given by 
each value of d b . Although these families of curves do not include all pos- 
sible shapes, one of them can usually be found to fit a practical shape for 
cushioning of this class over the anticipated range of use. 

Class D—Bi-linear Elasticity. This is characterized by a load-displace- 
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ment curve consisting of two straight line segments. The load displacement 
function is (see Fig. 1.4.4) 



P = k x 2 ^ .T 2 ^ d. 

P = hx 2 - {k b - k Q )d B x 2 ^ d a 



(1.4.4) 



It is useful especially in situations where very abrupt bottoming is possible. 

Class E — Hyperbolic Tangent Elasticity. When the mechanism of the 

cushioning is such as to limit the maximum force that can be transmitted 

over a considerable displacement range, the load-displacement function 



P = Pa tanh 



#0 #2 



(1.4.5) 



is useful. P is the asymptotic value of the force and k a is the initial spring 
rate (see Fig. 1.4.5). 





Fig. 1.4.5 

Fig. 1.4.4 — Bi-linear elasticity. Class D. 
Fig. 1.4.5 — Hyperbolic tangent elasticity. Class E. 

Class F — Anomalous Elasticity. In occasional instances the load-dis- 
placement curve of the cushioning cannot be matched accurately enough 
by any of the five preceding functions. In such cases a numerical integra- 
tion procedure can be used, as described in Section 1.15. 

1.5 Cushioning with Cubic Elasticity (Class B) 
Substituting (1.4.2) in (1.2.15) and performing the integration, we have: 



Now, let 



k ^J r r ^l = W , h . 
2 4 



/2W,h 



(1.5.1) 



(1.5.2) 
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that is, d is the displacement that would take place if the elasticity were 
linear (see equation (1.3.1)) with a constant spring rate k equal to the initial 
spring rate of the cubic elasticity. Also let 

4W 2 hr 



B = 



ki 



Then, from (1.5.1), (1.5.2) and (1.5.3) 



da 



-g(-i + Vi + b). 



:i.5.3: 



(1.5.4) 



Equation (1.5.4) is plotted in Fig. 1.5.1 which shows graphically how the 
maximum displacement d m compares with the "equivalent linear displace- 
ment Jo" as the parameter B is varied. Note that B depends on the weight 
of the packaged item, the height of drop and the shape of the load displace- 
ment curve (as determined by k and r). 




Fig. 1.4.6 — Anomalous elasticity. Class F. 

Similarly we can compare the maximum acceleration G m with the maxi- 
mum (Go) that would obtain if the load displacement curve were linear with 
spring rate k . The latter acceleration is given by 

Go = /§" d.5.5) 

and the former is obtained by finding P m from (1.2.17) and then, from 
(1.2.16), 

G m 



( y | (1 ^(^-ViV /; 



[1.5.6) 



Equation (1.5.6) is plotted in Fig. 1.5.2. 



1 .6 Procedure for Finding Maximum Acceleration and Displacement 
for Cushioning with Cubic Elasticity 

If the load-displacement curve of a cushioning system has the general 
appearance of Fig. 1.4.2 (where the slope increases or decreases gradually 
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with displacement) the following procedure may be used for estimating the 
effectiveness of the cushioning. 
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Fig. 1.4.7 — Family of load displacement curves for cushioning with tangent elasticity. 

a. Select the point on the load-displacement curve for which the load 
is equal to the weight of the packaged item multiplied by the allowable 
G m . Call this load P->. and the corresponding displacement d 2 . 
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Fig. 1.5.1 — Maximum displacement for cushioning with cubic elasticity. See 
equation (1.5.4). 
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Fig. 1.5.2 — Maximum acceleration for cushioning with cubic elasticity. See 
equation (1.5.6). 
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b. Select another point {d x , Pi) about half way toward the origin from 
(<£ 2 ,P 2 ). See Fig. 1.4.2. 

c. Calculate 

Pi ,2 P* .2 

7 02 " 7 o> 
fc = dl s J (1-6-D 



and 



di -di 



di d\ 
di — d\ 



(1.6.2) 



d. Using the known weight, W 2 , of the packaged item, the specified height 
of drop k, and k and r from (1.6.1) and (1.6.2), calculate B, d and Go from 
(1.5.2), (1.5.3) and (1.5.5). Then calculate the maximum acceleration G m 
and maximum displacement d m from (1.5.6) and (1.5.4) or find their values 
from Figs. 1.5.1 and 1.5.2. 

Example: A large vacuum tube, weighing 22.5 lbs, was packed in a 7" x 
7f" x 15" carton which was supported on corrugated cardboard spring pads 
in a 10!" x 11!" x 18 £" carton. The latter was, in turn, packed in 28 pounds 
of excelsior in a 25" x 25" x 30" carton. The tube is rated at 50g and the 
package is intended for a drop of three feet. 

A rod was inserted through a hole cut through the three cartons to the 
tube. Load was applied to the rod and the displacement of the tube was 
measured. The data obtained were 

P Xi 

(load in lbs) {displacement in inches) 

o 

100 fg 

200 I 

300 i 

400 IrV 

500 lft 

600 II 

700 11 

800 1ft 

900 If 

1000 1} 

1100 IB 

1200 13 

1300 1H 

1400 2 

The data are plotted in Fig. 1.6.1. The resulting curve is suitable for 
classification as either Class B or Class C cushioning. Considering it, for 
the present, as Class B, we take P 2 = 22.5 X 50 = 1225, and from the curve, 
di = 1.9 inches. Also, from the curve, take d\ = 1 inch and Pi = 365 
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lbs. Substituting these values in (1.6.1) and (1.6.2), we find k = 255 and 
r = 108. Then 

AWtkr 4 X 22.5 X 36 X 108 



B = 



K 



(255)' 



G ° = i/!°V 



2 X 36 X 255 



22.5 



= 5.4, 



= 28.6. 



Entering Fig. 1.5.2 with B = 5.4 we find G m /G = 1.9. Hence 

G m = 28.6 X 1.9 = 55 

This is close enough to the 50g rating of the tube to call the cushioning 
safe insofar as maximum acceleration is concerned. 
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Fig. 1.6.1 — Experimental load-displacement curve for a corrugated cardboard spring pad 

and excelsior cushion. 

The maximum displacement, obtained by entering Fig. 1.5.1 with B = 
5.4 and finding d m /d = 0.75. Then d m = 0.75 X 2h/G = 1.95 inches. 
Hence, the package is much larger than necessary since approximately 8 
inches of cushioning thickness is supplied to accommodate 2 inches of 
displacement. 

1.7 The Tension-Spring Package (Class B) 
The tension spring package is useful when the allowable G m is so small 
and height of drop so great that a large displacement (say d m > several 
inches) is required. The decision as to whether or not a tension spring 
package is indicated may be made on the basis of a preliminary estimate 
of displacement based on the linear case. Suppose the height of drop is 
to be 60 inches and the allowable acceleration for the packaged item is 
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20g. Then, from Equation 1.3.4, the approximate displacement that will 
be required is 

rf mS= 2X6 ° = 6 inches. (1.7.1) 

The actual maximum displacement in a tension spring package will prove 
to be somewhat more than 6 inches, but the preliminary calculation shows 
the displacement to be large enough to warrant the use of this type of 
cushioning. 

h g 




Q B 

Fig. 1.7.1 — Schematic diagram of a tension spring package. 

A schematic diagram of a typical tension spring package is shown in Fig. 
1.7.1 and a photograph of one design is given in Fig. 1.7.2. The packaged 
item is suspended on eight identical helical tension springs which diverge 
to the outer frame. The analysis and design procedures described in this 
and the following section apply equally well if the springs converge from the 
packaged item to the outer frame. With a slight modification, indicated 
in the next section, the procedure also applies if four of the springs (say, 
BJ, DL, EM, OG in Fig. 1.7.1) are omitted. 

In all cases, however, we shall consider only systems having reflected 
symmetry about each of three mutually perpendicular planes through the 
center of gravity of the packaged article. 
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Fig. 1.7.2 — A tension. spring package. 



The load-displacement characteristics of the spring system may be found 
by statical considerations. We shall examine, first, the displacement in the 
vertical direction in Fig. 1.7.1, using the following notations: 
P = force applied to the suspended object. 
.To = displacement of suspended object, 

.vo = perpendicular distance (IR, Fig. 1.7.1) from inner spring support 
point (/, Fig. 1.7.1) to nearest plane, perpendicular to displacement 
direction and containing four outer spring support points (A, B, C, 
D, Fig. 1.7.1); 
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li = distance (I A) between spring support points when suspended article 

is in equilibrium position, 
I = projection of l< on plane A BCD, 

f = I minus length (between hooks) of unstretched spring, 
k = spring rate of each spring. 

Consider, first, the action of one pair of springs, say EM and GO of Fig. 
1.7.1, independent of the remainder of the suspension. Since EM and GO 
lie in parallel vertical planes and the points M and remain in the initial 
planes of their respective springs during a vertical displacement, the two 
springs may be considered to lie in the same plane, and to be translated hori- 
zontally in this plane so that their outer ends are separated by a distance 21. 
Hence Fig. 1.7.3 may be used to represent the independent action of this 
pair of springs and it is required to find the force Q' needed to transform Fig. 
1.7.3(a) to Fig. 1.7.3(b). Initially there are two springs, each of length l — j 











V 
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Q 
^ a b 

Fig. 1.7.3 — Diagram used in discussion of tension spring package. 

and spring constant k, with no initial tension in them. One end of one spring 
is fixed at point E and one end of the other spring is fixed at a point G 
distant 21 from E. The springs are then stretched so that the two initially 
free ends are located at a point Y equidistant from E and G and distant x 2 
from line EG. The axis of each spring makes an angle a with EG, where 



sin a = 



x 2 



Vt 2 + 



*? 



In this state the axial force F in each spring is 



(1.7.2) 



(1.7.3) 



and the force Q', required to equilibrate the two forces F is IF sin a. Con- 
sidering the force Q' as a function of the displacement x% , we write 



2x',k 



<r ( 4) = -^=r 2 We + x?-t + f] 



(1.7.4) 



or 



where 
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1.7.5) 



'-?. 



»4 



Consider, next, the configuration shown in Fig. 1.7.4(a), where one end of 
each of four springs is fixed at a corner of a rectangle of length 21 and width 
2x . Each spring is again of length t — f. The four free ends of the springs 
are drawn together at a common point X at the center of the rectangle (see 
Fig. 1.7.4(b)). The system is in equilibrium in this position. A force Q 
is then applied at X in the plane of the rectangle and normal to the side 21. 



ejumjuuw> 



— / 






2X, 
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o/ X 2 



f f 

a b c 

Fig. 1.7.4 — Action of springs in a tension spring package. 

The common point X is displaced a distance x 2 to X' (see Fig. 1.7.4(c)). 
Writing z = x 2 /l, a = x /f, we observe that 

Q(z) = Q'(z + a) + Q'(z - a), (1.7.6) 

or, from equation (1.7.5), 

z + a z — 

(1.7.7) 



ew 



- 2M|2i-(l-*)["^= 



+ 



i-ayjj- 



+ ( 2 + aY ' Vl + (z 

The standard tension spring package has two sets of four springs so that 
the force P required to displace the common point X a distance x% is 

P(z) = 2Q(z). (1.7.8) 

If *2 is small in comparison with I (i.e., z is small in comparison with 
unity), equation (1.7.8) may be written approximately as 

(1- W 2-/)] (17p) 



P(z) = 4kC\ 2z - 



(1 + a 2 ) 3 ' 2 

Even when x 2 becomes almost as large as (, equation (1.7.9) has been found, 

experimentally, to be remarkably accurate. 

Writing 



^-(TT^J 



(1.7.10) 
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and 



U Tf 
C = 



(±±±L~ - i, (1.7.11) 

1 — 6 



equation (1.7.9) becomes 

p -«(« + 3s)- C 1 - 7 - 12 ) 

It is seen, by comparison with (1.4.2) that this is Class B cushioning 
(cubic elasticity). A' is the initial spring rate and c determines the rate of 
increase of stiffness with displacement. With the notation k , r oi Section 
1.5, we see that 

k = K (1.7.13) 

r - — . (1.7.14) 

Hence equations (1.5.6) and (1.5.4) may again be used to calculate maximum 
acceleration and displacement. B has the same meaning as before (Eq. 

1.5.3). 

To predict the performance, in the vertical direction (Fig. 1.7.1), of an 
existing tension spring package the same procedure as outlined in Section 
1.6 may be used, except that it is not necessary to have a load-displacement 
curve for calculating k and r. Instead, these parameters may be calculated 
directly from equations (1.7.10), (1.7.11), (1.7.13) and (1.7.14). The 
remainder of the procedure is the same as in Section 1.6(d). 

To predict the performance perpendicular to another face, say AEHD 
of Fig. 1.7.1, it is only necessary, in the calculation of k and r, to substitute 
.ro for .T , V for t (see Fig. 1.7.1) and, in place of b: 

b' - 1 - j, (1 - b). (1.7.15) 

The initial spring rate A" for any direction of acceleration may be calcu- 
lated from the initial spring rates A' x , A 2 , K 3 in the three directions normal 
to the faces of the frame by using the relation 

1 s 2 t- u- .. ,- .,. 

where 5, t, u are the direction cosines of the acceleration direction with 
respect to the normals to the faces of the frame. It is seen, from (1.7.16), 
that the spring rate is given by the radius to the surface of an ellipsoid whose 
principal semi-axes are A'i , A 2 , K 3 . 
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The displacement direction does not necessarily coincide with the acceler- 
ation direction. The angle 8 between them is given by 

-•-'K + K + fl (m7) 

where K is defined by equation (1.7.16). 

The spring characteristics may be made the same in all directions and the 
displacement direction may be made to coincide with the acceleration 
direction by setting 



and 

C= (' = I" 



V2 



(see Fig. 1.7.1). This makes b = b' = b", c = 0.828 and k/K = 0.274 
in the calculations of the next section. 

1.8 Procedure for Designing Tension Spring Packages 

The design of a tension spring package, as contrasted with the analysis of 
one, must proceed without initial knowledge of values for the parameters k 
and r, since these cannot be known until the springs are designed. There- 
fore equations (1.5.4) and (1.5.6) cannot be used directly. For design pur- 
poses they are transformed to the following set of formulas: 



v^^-^f-i 



(1.8.1) 



L = Vfe ~ I \/^K + VN + 8) 2 + f d.8.2) 

^J = V2(-l + Vf+5) (1-8.3) 

M 2 = N - ^ - | (1+ S)<-1 + y/lT+B) (1.8.4) 

7-*=-l + |/l+^. d.8.5) 

These formulas have been converted to design curves which are given in 
Figs. 1.8.1 to 1.8.5. The curves are for use in connection with the following 
routine procedure which has been found useful in designing the springs for 
tension spring packages. Reference should be made to Table I. 
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1. Enter, on Line 1, Table I, the weight (W 2 ) in pounds, of the sus- 
pended item. This includes the weight of the cradle or other holding 
arrangement and one-third the estimated weight of the springs. 



b=0.3 




0.1 0.2 0.3 0.4 0.5 6 0.7 0.8 9 1.0 



Xo 
X 

Fig. 1.8.1— Tension spring package design curve. Equation (1.8.1). 

2. Enter, on Line 2, the height of drop (h) in inches. 

3. Enter, on Line 3, the maximum allowable acceleration (G m ) in units 
of "number of times gravity." This should be determined before- 
hand from tests on the item to be packaged. 

4. Enter, on Line 4, the dimension x (inches). 
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5. Enter, on Line 5, the dimension I (inches). For a package to have 
the same spring rate in all directions, I = .ro\/2 is a necessary con- 
dition. 

6. Enter, on Line 6, the value chosen for b. As b becomes greater than 
zero, the stiffness of the whole suspension increases for a given stiff- 
ness of individual springs. The reverse happens for b less than zero. 



2.5 




Fig. 1.8.2 — Tension spring package design curve. Equation (1.8.2). 

7. Calculate xt>/t. 

8. Enter Fig. 1.8.1 with .To/^and find \/c. 

9. Calculate L = h/(^fc(G m ). 

10. Enter Fig. 1.8.2 with L and find N. 

11. Calculate K = {W£J)/(2hN). This is the initial spring rate of the 
suspension in the direction of x . 

12. Calculate/ = 3.13 (K/W 2 ) i . This is the natural frequency of vibra- 
tion (cycles per second) of the suspension for small amplitudes in the 
.To direction. This should not be close to the natural frequency of 
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vibration of any element in the packaged item, which should be 
determined by test beforehand (see, also, Section 4.2). In any case 
it is advisable to provide damping for the suspension. 



K 




Fig. 1.8.3— Tension spring package design curve. Equation (1.7.10). 

13. Enter Fig. 1.8.3 with x /£ and find k/K. If a four-spring package is 
desired, instead of an eight-spring package, (see Section 1.7) the 
value of k/K found on Fig. 1.8.3 should be multiplied by two before 
entering it on Line 13 in Table I. This is the only change required 
in the procedure. 

14. Calculate k = A'f - J. This is the spring rate of each of the springs in 
pounds per inch. 
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15. Calculate B = (2W,h)/(Kc( 2 ). 

16. Enter Fig. 1.8.4(a) or (b) with B and find d m /(Vct). 

17. Calculate djl - \/c-d m /y/c(. 

18. Calculated = l-djl. 

19. Calculate (d m /t) + (x /t). 
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0.4 0.6 0.8 10 1 ? 14 16 18 2.0 
B 
Fig. 1.8.4(a)— Tension spring package design curve. Equations (1.8.3) and (1.8.4). 

20. Enter Fig. 1.8.5 with {dj () + (x /C) and find (e/C) - b. e is the 
stretch of each spring (in inches) when the displacement is d m inches. 

21. Calculate F m = k-(c/ (')•(. This is the maximum load (in pounds) 
on each spring. 

22. 23, 24, 25, 26. These are the coil diameter, wire diameter, number of 
turns, fiber stress and length of coils. These quantities are calculated 
from the ordinary formulas, charts or slide rules for helical springs, 
using the values of k and F m from Lines 14 and 21. 

27. The length inside hooks is entered on Line 27 to group all of the spring 
specifications. 




0123456789 10 
B 
Fig. 1.8.4(b)— Tension spring package design curve. Equations (1.8.3) and (1.8.4). 
1.3 




Fig. 1.8.5— Tension spring package design curve. Equation (1.8.5). 
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t As an example of the calculations, Table I contains the entries for the 
design of springs for a 21 pound article (including | the estimated spring 
weight) which is to be packaged so as not to exceed 35g in a five-foot drop. 

TABLE I 

Computation Form for Tension Spring Packages 

1. W,(lbs.) 21 

2. h (ins.) 60 

3. G m 35 

4. .To (ins.) 5 

5. I (ins.) 7.07 

6. b 

7. Calc. xo/t 0.707 

8. Find \/c from Fig. 1.8.1 0.91 

9. Calc. h/y/c lG m = L 0.269 

10. Find N from Fig. 1.8.2 1 .265 

11. Calc. WiGm*/2hN = K (lbs/in.) 169.0 

12. Calc./ = 3.13 (A'/WolJ (cyc./sec.) 8.9 

13. Find k/K from Fig. 1.8.3 0.274 

14. Calc. k = K • k/K (lbs/in.) 46.5 

15. Calc. B = IWth/Kct 2 0.368 

16. Find d m /Vc I from Fig. 1.8.4 0.575 

17. Calc. djt = \/c ■ dJVcl 0.518 

18. Calc. d m = I • d m / 1 (ins.) 3.68 

19. Calc. d m /l + xo/t 1 .220 

20. Find e/t from Fig. 1.8.5 and line 6 0.580 

21. Calc. F m = k ■ e/ 1 • ^ (lbs) 191.0 

22. Coil diameter (ins.) 1 .40 

23. Wire diameter (ins.) 0.207 

24. Number of turns 19 

25. Fiber Stress (lbs./sq. in.) • 10 -3 80 

26. Length of Coils (ins.) 3.93 

27. Length inside books (ins.) 7.07 

1.9 Cushioning with Tangent Elasticity (Class C) 

This is one of the most frequently encountered classes of cushioning since 
it includes a very common type of bottoming (Figs. 1.4.3 and 1.4.7). The 
load-displacement function (equation (1.4.3)) takes into account the fact 
that the cushioning can be compressed only to a definite amount db . 

To find formulas for maximum acceleration and displacement, we pro- 
ceed as follows. Substitute equation (1.4.3) in (1.2.15) and perform the 
integration, obtaining 

4&o d b . ird m . . 

— — log cos — = -W 2 h, (1.9.1) 

■k l Idb 

which may be written as 
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Equation (1.9.2) can then be substituted into (1.4.3) to obtain the maximum 
force P m in accordance with (1.2.17): 



t-' — y^TBt)-*- ( ,3) 




Fig. 1.9.1— Curve for finding maximum acceleration for cushioning with tangent elasticity. 

See equation (1.9.4). 



The maximum acceleration is then obtained from (1.2.16) and may be 
written in the form 

Go ird 



j/^SJ 



- 1 



(1.9.4) 



where d and G are defined just as in (1.5.2) and (1.5.5). G is the maxi- 
mum acceleration that would obtain if the cushioning did not bottom, that 
is, if the spring rate remained constant at its initial value k . d is the 
maximum displacement that would be reached under the same linear con- 
ditions. Hence G m /G is a multiplying factor to be applied to a hypothetical 
linear cushioning to take into account the effect of bottoming. The multi- 
plying factor depends only on the ratio (d b /do) of the amount of space 
actually available to the amount of space that would be used under linear 
conditions. 
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The ratio G m /G is plotted against the ratio d b /d in Fig. 1.9.1. It may- 
be seen that the multiplying factor increases very rapidly as the displace- 
ment ratio (db/do) falls below unity. For example, if the cushioning, with 
tangent elasticity, reaches hard bottoming (d b ) when only 80% of the 
required displacement (d ) is attained, the acceleration is multiplied by 
3.5; if only 60% of the required displacement is available, the acceleration 
is multiplied by 11.5. 

Example: To illustrate with a numerical example, consider the case already 
discussed in Section 1.3, where we found that a spring rate of 694 lbs/in 
and a displacement of 1.44 inches were required to limit a 20-pound article 




.2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 24 2.6 2.8 3.0 



Fig. 1 .9.2 — Curve for finding maximum displacement for cushioning with tangent elasticity. 

See equation (1.9.5). 

to an acceleration of 50g in a 3-foot drop with linear cushioning. Let us 
suppose that only 1.15 inches are available, instead of 1.44 inches, and that 
the cushioning has tangent elasticity starting with a spring rate of 694 
lbs./in. Entering the curve of Fig. 1.9.1 at d b /d = 1.15/1.44 we find 
G„,/Go = 3.5. Hence the maximum acceleration will be 175g instead of the 
required 50g. This illustrates the wide variations in acceleration that may 
occur as a result of minor variations of dimensions in high G packages. 

It is not necessarily true that the 175g test is 3.5 times as severe as the 
50g test for all elements of the supported structure, since the severity de- 
pends also on the shape and scale of the acceleration-time relation. The 
factor may be more or less than 3.5 but it will be very close to this value for 
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all high-frequency elements of the structure. This subject is treated in 
detail in Parts II and III. 

The maximum displacement d m , in the case of tangent elasticity, may be 
calculated from equation (1.9.2) or, in terms of d b /d , from 



d m 2 _i 

— = - cos 
di, v 



-[-T®1 (1 -°- 3) 



The ratio djd b is plotted against d b /d in Fig. 1.9.2. 

The use of Fig. 1.9.2 can be illustrated with the example already calcu- 
lated, in which d b /d u = 1.15/1.44 = 0.8. Entering the abscissa of Fig. 
1.9.2 with d b /do = 0.8 we find d m /d b = 0.915. Hence the maximum dis- 
placement will be 0.915 X 1.15 = 1.05 inches. 

1.10 Optimum Shape of Load-Displacement Curve for Tangent 

Elasticity 

It is possible to choose the best shape for the load-displacement curve 
of the cushioning from those represented in Fig. 1.4.7. This will be, of 
course, not the best of all possible curves, but only the best among "tangent 
elasticity" curves. The best shape is defined as the one that yields the 
smallest maximum acceleration (G m ) for a given weight (IF 2 ), height of drop 
(h) and available space (d b ). This leaves the initial spring rate (* ) as the 
only remaining variable. To find its optimum value (say k ), set equal to 
zero the derivative of G m (equation (1.9.4)) with respect to k , remembering 
that Go and d are functions of k . The result is 

/ i*W*h \ /VIF 2 /A 

O-l^V^tooj-^ ' (U01) 

from which 

d b 

Substituting (1.10.2) in (1.9.4) we find the minimum value (G' m ) of maxi- 
mum acceleration to be 

G' m = 3 -^. (1-10.3) 

d b 

To illustrate the application of equations (1.10.2) and (1.10.3), consider 
again the case of the 20-pound article dropped from a height of three feet. 
We found that a linear spring, with a spring constant of 694 lbs/in, would 
limit the maximum acceleration to 50g if 1.44 inches of displacement were 
available. If only 1.15 inches of displacement are available, and the initial 
spring rate is kept at 694 lbs/in, we found the maximum acceleration to be 
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175g if the cushioning bottoms with tangent elasticity. Now, according to 
equation (1.10.2), the best initial spring rate for cushioning with tangent 
elasticity would be 

,, 3.1 X 20 X 36 .,__., ,. 
= (1 151 2 = U»/in. 

In this case, equation (1.10.3) gives, for the maximum acceleration, 

, 3.9 X 36 

Gm = "Tis" = 122g ' 

Hence, confronted with a space limitation less than that required for a 50g 
linear spring, it is better to use an initial spring rate higher than that for the 
50g linear spring in order to strike an economical balance between displace- 
ment and bottoming. The best balance, among cushionings having tangent 
elasticity, is obtained by using equation (1.10.2). 

If no factor of safety is considered, it would be still better not to use a 
bottoming type of cushion at all. From equations (1.3.5) and (1.3.3) 
it can be seen that a linear spring with a constant of 1090 lbs/in will give 
only 63g with a displacement of 1.15 inches. Such a spring, though, would 
bottom very sharply at a drop slightly higher than 3 ft. and would give 
an acceleration much greater than cushioning with tangent elasticity which 
bottoms more gradually. This may be important if there are high-fre- 
quency, brittle elements in the packaged article (see Part III). 

1.11 Procedure for Finding Maximum Acceleration and 

Displacement for Cushioning wtth Tangent Elasticity 

(Class C) 

To illustrate the use of the equations and curves for Class C cushioning, 
the same example used for Class B will be used, as it was observed that the 
experimental load-displacement curve in that example (Fig. 1.6.1) is a 
border line one which can be treated as either B or C. 

By laying a straight edge along the first part of the curve (Fig. 1.6.1), 
the average initial spring rate is found to be 305 lbs/in. This value is taken 
as k in the present case. 

The next step is to find a value of db such that a graph of P/db vs x^/db 
will fall slightly above the curve k = 30(0) in Fig. 1.4.7; d b must be greater 
than 2 inches, since that displacement was obtained in the experiment. 
As a trial take d b = 2.25 inches and test it at one point, say the experi- 
mental point P = 300 lbs., .v 2 = \ in. Then P/db = 133 and x 2 /db = 
0.39. The point (0.39, 133) falls below the curve k = 30(0) in Fig. 1.4.7. 
Next try db = 2.5 inches. In this case, for the experimental point P = 
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300, x 2 = |, we find P/d b = 120, x 2 /d b = .35. This falls slightly above the 
k = 30(0) curve as required. The whole experimental curve is then 
plotted to the coordinates P/2.5 vs Xj/2.5 and is found to fit as closely as 
necessary. Hence the parameters are adopted as k = 305, d b = 2.5. 

We can now calculate the maximum acceleration that the tube will 
receive in, say, a three-foot drop test. First calculate, from equations (1) 
and (2), 



-^ = \/ 



2 X 36 X 22.5 = „ -j 
305 



Go m ym ' y 



2 X 36 X 305 = 31 3 



Then d b /d = 1.08. Entering Fig. 4 with this value we find G m /G = 1.82. 
Hence the maximum acceleration is: 

G m = 31.3 X 1.82 = 57g. 

Finally, entering Fig. 5 with d b /d = 1.08 we find d m /d b = 0.8. Hence the 
maximum displacement is d m = 0.8 X 2.5 = 2.0 inches. This indicates 
that the load-displacement test was carried far enough to cover the range 
up to a three-foot drop. 

It may be observed that the results obtained, by treating the same data 
as Class B or Class C cushioning, agree within a few per cent. This is 
because, in the example chosen, both B and C curves can be made to fit the 
experimental load-displacement curve. 

1.12 Consequences of Abrupt Bottoming (Class D) 
It is useful to examine cushioning systems that can bottom more abruptly 
than Class C cushioning. Abrupt bottoming is possible, for example, in a 
tension spring package lacking a snubbing device. An estimate of the 
increase in acceleration can be made by studying the case of bilinear elasticity 
(Fig. 1.4.4). Here we have a spring rate k up to a displacement d„ follow- 
ing which the cushioning has a different spring rate k b . k represents the 
average spring rate before bottoming and k b can represent the much greater 
stiffness of the wall of the container. 
If d > d s , that is, if 

^>d 8 , (1.12.1) 

the suspended article will bottom and the maximum displacement and 
acceleration are obtained by using both of the equations (1.4.4) in evaluating 
the integral in (1.2.15). Thus, 

'*' k x, dx 2 + /"'" [to - (h - h)d s ] dx, = W,h. (1.12.2) 



i 



[ 

Jo 
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The remainder of the procedure for finding G m is the same as before. 
The value of d m found from (1.12.2) is substituted for ;v 2 in the second of 




Fig. 1.12.1 — Curves for finding maximum acceleration as a result of abrupt bottomine 

See equation (1.12.3). °' 

(1.4.4) and the value of P m , thus obtained, is used, in (1.2.16), with the 
result: 

where G is the acceleration that would be reached if a displacement d 
were available: 



r /2hko 



(1.12.4) 



The ratio G m /G is plotted against d s /d in Fig. 1.12.1 for several values 
of k b /k . Since, in practice, k b might be thousands of times as great as k Q , 
it may be seen that the increase in maximum acceleration can be very large 
even when d a is only slightly less than d . It is apparent that a snubbing 
device is desirable in a tension spring suspension. This is especially true 
when considering high-frequency elements of the packaged article. It will 
be shown, in Part III, that low-frequency elements are not affected as much 
as might be expected from consideration of maximum acceleration alone. 
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1.13 Cushioning with Hyperbolic Tangent Elasticity (Class E) 

In the preceding sections, there have been considered four types of elas- 
ticity (linear, cubic, tangent and bilinear) that fit the load-displacement 
characteristics of the more common cushioning materials and devices. 
There now remains the problem of finding more nearly ideal shapes of 
elasticity. By "more nearly ideal" is meant a shape which will result in a 
smaller maximum displacement for a given maximum acceleration. This is 
important in the packaging of very delicate articles if shipping space is 

limited. 

It may be observed (equation (1.2.15)) that the total area under the 
load-displacement curve is equal to the maximum energy of the system. 
The maximum ordinate of the enclosed area is proportional to the maxi- 
mum acceleration. Hence, if we wish to (1) limit the maximum acceleration 
(2) accomodate a given kinetic energy and (3) have as small a displace- 
ment as possible, the best shape for the load displacement function is P = 
constant, where the constant is the product of the supported mass and the 
maximum allowable acceleration. 

It is not practical to obtain this ideal shape exactly, for there will always 
be a finite initial spring rate and a rounding off of the load-displacement 
curve to the limiting maximum load. A function which represents this 
practical condition (and also includes the ideal case) is the hyperbolic 
tangent function mentioned in Section 1.4: 

P = P tanh ^p- (1.13.1) 

The formulas for maximum acceleration and displacement are found in 
the same way as for the other classes of cushioning with the results: 

<C= 7 ;°cosh-exp(^) (1.13.2) 



or 



and 



or 



,/ .^cosh-'exp^^l) (1-13.3) 



Wt Go 



Po, , h d m 



G m =^-°tanh^ (1.13.4) 



a.-Sy^-q?) 
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where, as before 



*-V?, *-</ 



2hk 



Equations (1.13.3) and (1.13.5) are plotted, in Figs. 1.13.1 and 1.13.2, 
against the dimensionless parameter P /W 2 G . The latter is the ratio of 
the maximum force, that the hyperbolic tangent cushioning will transmit, 




WoGr 



Fig. 1.13.1 — Maximum displacement for cushioning with hyperbolic tangent elasticity. 

See equation (1.13.3). 

to the force that linear cushioning would transmit under the conditions 
specified. 

To find the value of ko which yields the minimum value of acceleration 
for a given maximum displacement, differentiate (1.13.4) with respect to k 
and set the result equal to zero: 



sech 2 ^-" = 0. 



(1.13.6) 



This is satisfied by k — > » , which represents the rectangular load dis- 
placement curve and confirms the conclusion reached from energy 
considerations. 
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Taking the limit of (1.13.4) as fo -» oo , we find the optimum acceleration 
to be 






(1.13.7) 



The corresponding maximum displacement is found, from (1.13.2) to be 

£«^-£- d-13.8) 



P Gm 




W 2 G 



Fig. 1.13.2— Maximum acceleration for cushioning with hyperbolic tangent elasticity. 

See equation (1.13.5). 

1.14 Minimum Space Requirements for Various Classes or 

Cushioning 

It is interesting to compare the minimum amount of space for displace- 
ment that can be attained with the various kinds of cushioning that have 
been discussed. 



Hyperbolic Tangent Elasticity 
Linear Elasticity 
Tangent Elasticity 



dm 


= 


h 

2h 


d m 


= 


G m 



dm = 



3.9A 
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Cubic elasticity will give a d m somewhat more or less than 2h/G m depending 
upon whether the parameter r is positive or negative. 

It is seen that a factor of almost four can be gained, in the linear dimensions 
of the cushioning space required, by replacing the tangent type of cushioning 
with the hyperbolic tangent type. 

There are several ways of obtaining a load-displacement curve with a 
shape similar to the hyperbolic tangent curve. One of the most interesting 
is suggested by the fact that the load-displacement curve of a strut has 
approximately this shape. Hence a bristle brush has the proper 
characteristics. 
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0.100 
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49.8 
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9.2 


31.6 


7 


0.050 


1.050 
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32.9 
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11.1 


39.5 
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0.050 
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950 


42.0 


245.8 


13.3 


51.4 
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0.050 
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1370 
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74.0 
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0.025 
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91.0 


11 


0.025 


1.200 


2240 


49.0 


390.9 


21.1 


121.0 


12 


0.0125 


1.2125 


2620 


30.4 


421.3 


22.8 


141.5 


13 


0.0125 


1.225 


3200 


36.4 


457.7 


24.7 


173.0 



1.15 Numerical Method for Analyzing Class F Cushioning 
The numerical method to be described is one that has been adapted from 
a graphical one used by the Committee on Packing and Handling of Radio 
Valves of the British Radio Board. The method has advantages of sim- 
plicity in concept and ease of application, especially when the load-displace- 
ment curve of the cushioning does not resemble closely one of the Classes A 
to E described above. It has the disadvantage that it does not yield, 
directly, numerical factors by which the spring rate or depth of cushioning 
should be changed in the event that the analysis reveals inadequate or more 
than adequate protection. 

The method is based on the fact that the area under the load-displace- 
ment curve of the cushioning represents the energy stored in, or absorbed 
by, the cushion. The total amount of energy that must be transferred is 
equal to the product of the weight (W 2 ) of the suspended item and the 
height (//) of drop. By finding the abscissa (.v 2 ) and its ordinate (P) which 
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include an area W 2 h, the maximum displacement is immediately x 2 and the 
maximum acceleration is the quotient P/W 2 , in accordance with equations 
(1.2.15) and (1.2.16). 

As actually applied in the present instance, the British method was 
modified slightly to make the procedure a routine numerical one. The 




0.1 0.2 0.3 



0.4 0.5 0.6 0.7 0.8 0.9 
DISPLACEMENT (INCHES) 



1.0 1.1 1.2 



Fig. 1.15.1 — Experimental load-displacement curve for Table II. 

computing form is given in detail in Table II, in which the data are taken 
from an experimental load-displacement curve (Fig. 1.15.1) for the end spring 
pads of a vacuum tube package. The load-displacement data are listed in 
Columns 3 and 4 of Table II. The meaning of each column in the table is as 

follows. 

Column 1. n(= 1, 2, 3- • •) is the number that identifies the displacement 
(and corresponding load) up to which the area under the load- 
displacement curve is to be calculated. 
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Column 2. A(x 2 ) n is the increment of displacement between (* 2 )„_i and 
(x 2 ) n . A(> 2 )„ = (*»)■ - Wn-i , see Fig. 1.15.2. Note that, 
as the curve becomes steeper, A(.*r 2 )„ is taken smaller for better 
accuracy. 

Column 3. (.r 2 )„ is the displacement associated with the n point (see Fig. 

1.15.2). 
Column 4. P„ is the load that produces displacement (.r 2 ) n . 




A(X 2 ) t 



(X 2 ) 2 (X 2 ) 3 (X 2 ) 4 ,\p<2) n 

*+* *\+ — *| H-H~< x 2Xm 

V A(X 2 ) 



A(X 2 ) 2 A(X 2 ) 3 ATX 2 J 

X 2 = DISPLACEMENT 



Fig. 1.15.2 — Graphical illustration of numerical method of calculating area under load- 
displacement curve. See Table II. 



Column 5. AA n = %A(x 2 )»(P n -i + P n ) is the area of the trapezoid with 
altitude A(.r 2 )„ and bases P„_i and P„ . It is approximately the 
energy absorbed by the cushioning in displacing from (# 2 ) n _i 

to (x 2 ) n ■ 

Column 6. A n is the sum of all the trapezoidal areas from x 2 = to x 2 = 
(x 2 )„ . It is approximately the total energy the cushioning can 
absorb in displacing an amount (x 2 ) n beginning at zero dis- 
placement. Note that A is always equal to zero. 

Column 7. //„ = A„/W 2 is the height of fall that will cause the cushioning 
to displace an amount (x 2 ) n . In Table II, W 2 = 18.5 pounds. 

Column 8. G„ = P n /W 2 is the maximum acceleration experienced by the 
suspended mass when dropped from a height //„ . 
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HEIGHT OF DROP (INCHES) 

Fig. 1.15.3 — Maximum acceleration vs. height of drop for an 18.5 pound article supported 
on cushioning with the load-displacement curve of Fig. 1.15.1. Sec Table II. 

From the last two columns of the table a curve of height of drop vs. the 
corresponding acceleration may be plotted as in Fig. 1.15.3. 

PART II 
ACCELERATION-TIME RELATIONS 

2.1 Introduction 

In Part I we were concerned primarily with the maximum acceleration of 
the packaged item. In this part we shall study the details of the variation 
of acceleration with time in order to have this information available for our 
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study, in Part III, of its influence on the response of elements of the packaged 
item. 

The first case to be considered will be the simple single mass and linear 
spring example described in Sections 1.2 and 1.3. Following this the 
phenomenon of rebound of the package will be considered. The influence 
of velocity damping and dry friction will be studied; and, finally, the effects 
of non-linearity of the cushion elasticity on the acceleration-time relation 
will be investigated. 

2.2 Acceleration-time Relation for Linear Elasticity 

Returning to the elementary example studied in Sections 1.2 and 1.3, we 
first write the equation of motion for the mass m 2 , on its linear spring of 
spring rate k 2 (see Fig. 1.2.1.). Equation (1.2.3) becomes 

m 2 x 2 + k 2 x 2 = m-ig. (2.2.1) 

Using the initial conditions 

Nm = 0, (2.2.2) 

[x 2 ]t=o = y/lgli, (2.2.3) 

the solution of (2.2.1) is 



Vg 2 + 2w;gA . , , N . g 

= — - — 4 — sin (o> 2 / — a) + -2 



* 2 = -^-s — ^ — — sin (« 2 / — a) + -^ (2.2.4) 



or 



t /m 



4- d- m sin (a) 2 / - a) + I- 2 , (2.2.5) 

«2 



where 



and 



- - l/l - ** - 1 < 2 - 2 - 6 > 



e Wo 

a = tan" 1 7= = tan" 1 — (2.2.7) 

to 2 V2gh k 2 d m 

wo is the circular frequency, f 2 is the frequency and T 2 is the period of vibra- 
tion of the mass m 2 on its spring; d„, has the same definition as in Section 1.3 
(equation (1.3.1)). 

Now, W 2 /k 2 is the static displacement of the mass m 2 on its spring. This 
is usually very small in comparison with the maximum displacement (d m ) 
during impact. Hence H'2/^2 will be neglected, and (2.2.5) becomes 

X» = d m sinw 2 t. (2.2.8) 
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Differentiating (2.2.8) twice with respect to t, we find, for the acceleration 
x 2 = —uldm sin u 2 t = — m s/2gh sin w 2 /. (2.2.9) 

Hence the absolute magnitude of the maximum acceleration is 

r | Xj | max _ &2 d m _ / 2hki ,« 9 10) 

T~ g ~ V Wi K " * 

as before. 





Fig. 2.2.1 Fig. 2.2.2 

Fig. 2.2.1— Half -sine-wave pulse acceleration. See equation (2.2.9). 
Fig. 2.2.2— Oscillogram of a half-sine-wave pulse obtained with a piezo-crystal 

accelerometer. 

Equation (2.2.9) shows that the acceleration varies sinusoidally with time. 
It rises from its initial zero value to its maximum in a time 7r/2co 2 , at which 
time the displacement also reaches its maximum value. The acceleration 
returns to zero again at time ir/u 2 . At this time the displacement is also 
zero. This is the end of the range of applicability of equation (2.2.9); for 
when I becomes slightly greater than 7r/a> 2 , a tension in the spring is required. 
Since no mechanism, such as a large mass m 3 (Fig. 0.2.2), has been supplied, 
to allow a tension in the spring to develop, the system will rebound from the 
floor at the end of the half period 7r/w 2 . The acceleration is thus a half- 
sinusoidal pulse of duration ir/o> 2 = T 2 /2 and amplitude G m g as illustrated 
in Fig. 2.2.1. An oscillogram of such a pulse obtained with a piezo-crystal 
accelerometer is shown in Fig. 2.2.2. 

2.3 Package Rebound. 
The presence of the mass of an outer container will affect the acceleration 
after the first half cycle of displacement. The outer container is represented 
by the mass Wj in the general idealized system illustrated in Fig. 0.2.2 and in 
the simpler system (Fig. 2.3.1) that we shall consider now. 
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Two pairs of equations are necessary to describe the action of the system; 
one pair applies during the time of contact of m 3 with the floor and the 
second pair applies if rebound occurs. 

The mass m 3 is assumed to be inelastic (see Section 0.2) so that, during the 
interval of its contact with the floor, the equation of motion for m 2 will be 
the same as before (2.2.1). In addition, there will be an equation of equili- 
brium for the mass m 3 : 



R = k 2 x 2 + m 3 g , 
where R is the upward force exerted by the floor on m 3 . 



(2.3.1) 




h 



m 3 



77777777T77777 



Fig. 2.3.1— Two-mass system representing packaged wticle, linear cushioning and 

outer container. 

Equations (2.3.1) and (2.2.1) will hold as long as .R is positive. To find out 
when R > 0, solve (2.2.1) for k 2 x 2 and substitute in (2.3.1): 



R = W 2 + W 3 - m 2 x 2 . 



(2.3.2) 



That is, a necessary condition for rebound is that the mass of the cushioned 
article, multiplied by its maximum acceleration, exceeds the total weight 
of the package. The condition for rebound may be written 



G m > 



W 2 + W 3 
W 2 



(2.3.3) 



This is a necessary, but not a sufficient, condition for rebound because there 
will be energy losses as a result of damping and permanent deformation. 
G m will generally have to be considerably greater than the right hand side 
of (2.3.3) for rebound to occur. 

If rebound does not occur, equation (2.2.9) continues to apply, except for 
damping which will be considered in Section 2.5. 
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2.4 Motion After Rebound 
If rebound occurs, the equations of motion for the two masses, tit* and m 3 , 
are 

»*2.i : 2 + ki(x2 — x 3 ) = tthg, (2.4.1) 

m 3 x 3 - k 2 (xi - x 3 ) = Wag. (2.4.2) 

Multiplying (2.4.1) by m 3 and (2.4.2) by m 2 and subtracting, we find 

my + k 2 y = 0, (2.4.3) 

where 

y = x% — x 3 , 

mim 3 



m = 



W2 + *»3 



(2.4.4) 
(2.4.5) 



^^xy ^ 



r\ a /> « ** *• *■ 



Fig 2 4 1— Oscillogram illustrating the half-sine pulse followed by the higher frequency, 
lower amplitude vibration of the packaged article in a rebounding package. 

Equation (2.4.3) is the equation governing the vibration of the two-mass 
system as a simple oscillator. The circular frequency of the vibration is 

co = A A 2 (2-4.6) 

y m 

and it may be noticed that this frequency is always greater thanco 2 (equation 
(2.2.6)). This fact is important in estimating the effect of vibrations on 
elements of the packaged item (Section 3.5). 

co is also the frequency of vibration of the packaged article during the 
interval of free fall. This vibration (usually of small amplitude) is initiated 
by the sudden release of the dead load displacement of the packaged article. 

As an intermediate step in obtaining the acceleration after rebound we 
shall find the magnitude of the relative displacement (y) of the two masses. 
To do this it is necessary to solve equation (2.4.3) with the appropriate 
boundary conditions. Calling U the time at which m 3 leaves the floor, we 
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must find y and y at / = t T . Since m 3 is motionless at t = t r , the relative 
displacement and velocity at that time are identical with x 2 and .r 2 respec- 
tively. The former is simply the stretch of the spring necessary to just 
pull the mass m 3 off the floor, i.e., 

M« r = [**U r - ~. (2.4.7) 

To find the velocity at / = t T , substitute (2.4.7) in (2.2.4) and also substi- 
tute t r for t in the latter. This gives an equation for determining t r . Then, 
returning to (2.2.4), differentiate it once to obtain x 2 and substitute for I the 
value t T just found. The result is 



[*Ji-« r = [y]t=t r = - A/ 2gh - 



gir^n'o + ny^ {2AS) 



The solution of (2.4.3) with initial conditions (2.4.7) and (2.4.8) is 



y = - ~ A/2gh -^f 8 sin («/ - f ), (2.4.9) 

where f = co/ r - tan" 1 #^ . 

[y]t=t r 

We are now in a position to find the acceleration of the packaged item 
after rebound. Substitute y of (2.4.9) for * 2 — * 3 in (2.4.1) to obtain 

£2 = g+- A/lgh - T %- g sin (cot - f). (2.4.10) 

co y ki 

To obtain a simple formula for the ratio of the maximum accelerations 
after and before rebound, let us assume that both are much greater than 
gravitational acceleration. Then if 

G r = maximum number of g's after rebound, (maximum of (2.4.10)) 
Gm = A/ -ttt = maximum number of g's before rebound, 
we find, from (2.4.10), neglecting the term g outside the radical, 



Gr ' l! /l-^. (2.4.11; 



G m W 2 +W*y 2hh 

Hence, the maximum acceleration after rebound is always less than the 
maximum acceleration before rebound. Therefore, conditions after re- 
bound need only be examined when the frequency after rebound (see equa- 
tion (2.4.6)) is near the natural frequency of vibration of a critical element 
of the packaged item (see Section 3.5). 
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The complete acceleration history of a rebounding package with un- 
damped linear cushioning is thus a half sine wave pulse of amplitude G m = 
y/lhkijWi and duration t/u 2 followed by an oscillating acceleration of 
amplitude given by (2.4.11) and frequency given by (2.4.6). Such a wave 
shape is shown in Fig. 2.4.1. 

2.5 Influence of Damping on Acceleration 

The presence of damping in cushioning is always desirable to prevent the 
building up of large amplitudes as a result of periodic disturbances. How- 
ever, damping also has an effect on the maximum acceleration that is at- 
tained in a drop test. From the latter point of view there is an optimum 
amount of damping and an amount that should not be exceeded if the maxi- 
mum undamped acceleration is not to be exceeded. 

We shall consider the case of a linear cushion with damping proportional 
to velocity. The system is represented in Fig. 2.5.1. With the addition 



— //////////sT?/ // ' ' 
Fig. 2.5.1 — Idealization of linear cushioning with velocity damping. 

of the damping term the equation of motion of m* , during contact of the 
package with the floor, is 

W2»2 + ^2*2 "T" ^2*2 = 0, (2.5.1) 

in which c% is the damping coefficient of the cushioning. Equation (2.5.1) 
is more conveniently expressed as 

Xi + 2/32C02-i"2 + &>2#2 = 0, (2.5.2) 



where 



0)2 = 



TV 



(2.5.3) 



2W20>2 

w 2 is the undamped circular frequency of vibration of m 2 on its spring and /3 2 
is the fraction of critical damping. /3 2 = means no damping and /S 2 = 1 
means just enough damping so that there will be no oscillation if the pack- 
aged article is displaced and released. 
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The acceleration solution of (2.5.2), with the initial conditions of the drop 
test (see (2.2.2) and (2.2.3)) is 



where 



_ _ <aiV2gh 0*- 

•1-2 — — / ^. t. 

vi - /?! 



tan 7 = 



cos (« 2 <Vl - PI + y) (2.5.5) 

2/3| - 1 



2/Wl ~/3|* 



(2.5.6) 
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Fig. 2.5.2 — Acceleration-time curves for linear cushioning with various amounts of 
damping (no rebound). See equation (2.5.5). 

The acceleration is thus a damped sinusoid with an abruptly reached 
initial value whose magnitude depends upon the amount of damping. For 
small damping, the initial acceleration is small and then the acceleration 
increases, but never reaches the value that would be reached without any 
damping. For high damping (fit > 0.5) the initial value is greater than 
without any damping and falls off thereafter. Figure 2.5.2 shows the shapes 
of the acceleration time curves for several values of /3 2 . All of the curves 
are for no rebound. It may be seen, from equation (2.5.5) and Fig. 2.5.2 
that the addition of damping changes the shape of the acceleration-time 
relation in three ways. First, a damped sinusoid replaces the pure sinusoid; 
second, the frequency is reduced; and, third, the initial phase is changed. 

It is useful to consider in detail the effect of damping on maximum accel- 
eration. Let 

G m = maximum number of g's with damping 

Go = A/ t=t- 2 = maximum number of g's without damping. 

V w 2 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
fi 2 (FRACTION OF CRITICAL DAMPING) 

Fig. 2.5.3— Influence of velocity damping on maximum acceleration. 



Then, from (2.5.5), at / = 



and, after / = 0, 



Go 



= 2fa 



See equation (2.5.5). 



(2.5.7) 



GL. 

G =e 



_ a -Pvtt, 



(2.5.8) 



where /„, , the time at which the maximum occurs, is given by 

tan a> 2 /„, Vl - A ^ 2 (3 - 4j 8l) " ( j 

The largest value of G m /G from (2.5.7) and (2.5.8) is plotted against fa 
in Fig. 2.5.3. It is shown there that, as the damping is increased from zero, 
the maximum acceleration first decreases to a minimum of 80% of G and 
then increases to G at 50% of critical damping. In this interval the maxi- 
mum acceleration occurs after t = 0. For damping greater than fa = 
0.5 the maximum acceleration occurs at the instant of contact and increases 
in direct proportion to fa . 

2.6 Influence of Damping on Rebound 

In considering rebound without damping, it was found that rebound does 
not occur unless the product of the maximum acceleration and the sus- 
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pended mass exceeds the total weight of the package. It was not necessary 
to distinguish between maximum acceleration on the first downstroke and 
first upstroke, since these are the same when there is no damping. With 
damping, however, the maximum acceleration on the first downstroke is 



*§ 

















































































1.0 


























\ 
























0.8 


\ 
























\ 
























0.6 


\ 


V 
























\ 






















0.4 




\ 














































0.2 













































































1.0 



fli 



Fig. 2.6.1 — Influence of velocity damping on maximum upstroke acceleration. See 

equation (2.5.5). 



greater than that on the first upstroke (Fig. 2.5.2) and it is the latter that 
controls rebound. Hence damping inhibits rebound. 

For example, with 50% of critical damping (/3 2 = 0.5), equations (2.5.8) 
and (2.5.9) and Fig. 2.5.2 show that for the first downstroke G,„/Go = 1 
while for the first upstroke G m /G = 0.164. Hence the tendency to rebound 
is reduced by a factor of six when damping to the extent of 50% of critical 
is added to an undamped package. 

The ratio of the maximum acceleration on the first upstroke to the maxi- 
mum undamped acceleration is plotted in Fig. 2.6.1 for various values of /3 2 . 
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2.7 Influence of Dry Friction on Acceleration and Displacement 
By "dry friction" is meant friction that is independent of velocity except 
for sign. During contact of the package with the floor the motion of m 2 
might be opposed by a constant friction force F. Such a force is developed, 
for example, in a package with corrugated spring pad cushioning by rubbing 
against the side and end pads in a top or bottom drop. A typical idealized 



2F 




Fig. 2.7.1— Load vs. displacement for cushioning with dry friction. 

load-displacement curve is shown in Fig. 2.7.1. For the first downstroke 
of m 2 , the equation of motion of m 2 is 

m 2 x 2 + k 2 x 2 = —F. (2.7.1) 

With initial conditions 

[* 2 ] (=0 = 0, fo]«i = y/2gh (2.7.2) 

the solution of (2.7.1) is 



Xa 



where 






(2.7.3) 



<^/ 2j P> 



tan a = 



F 



k 2 d W 2 Go ' 
n . /2hk 2 

Go - y w 2 ■ 



DYNAMICS OF PACKAGE CUSHIONING 405 

d and Go are the maximum displacement and acceleration that would obtain 
if no friction were present. From (2.7.2) the maximum displacement with 
friction is 

*--|/« + (3'-5- (2 - 7 - 4) 

Hence, the presence of friction decreases the maximum displacement since 
d m < do. 

From (2.7.3) the acceleration is 



= " jAo 2 + (j0 sin («*' + «). (2-7.5) 



so that the maximum acceleration is 

G- = j/cS + (0, (2.7.6) 

which is greater than the maximum acceleration without friction. 

It would appear, at first glance, that cushioning with friction always 
gives a greater acceleration than the corresponding cushioning without 
friction. However, the reverse is actually true provided we allow the same 
displacement in both cases. This may be done, as may be seen from (2.7.4), 
by decreasing the spring rate in the cushioning with friction to 

2F 
k, = kt-=-. (2.7.7) 

Oo 

The maximum acceleration in the cushioning with friction is then, from 
(2.7.6), 

&-<*-£, «*w*% (2 - 7 - 8) 

That is, for the same maximum displacement, the maximum acceleration 
is reduced by the addition of dry friction. 

2.8 Acceleration-Time Relation for Cubic Elasticity 

As an example of the effect of nonlinearity of the cushioning on the shape 
of the acceleration-time function, the case of cubic elasticity (Class B) 
will be analyzed. The system to be considered is illustrated in Fig. 1.2.1, 
and the load-displacement relation for the cushioning is given by 

P = kox 2 + rxl . (2.8.1) 
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Substituting (2.8.1) in (1.2.13) and performing the indicated integration, 
we find 



.2 r, , k 2 r 4 

X 2 = 2gk — .To — — x 2 . 

Remembering that x 2 = dx 2 /dt, we solve (2.8.2) for dt: 

dx 2 



dt = 



• 



2^-^-x\ -/-* 2 



(2.8.2) 



(2.8.3) 
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Fig. 2.8.1— Duration of acceleration pulse for cushioning with cubic elasticity. See 

equation (2.8.10). 



Then, with the initial condition ;v 2 = when / = 0, the integral of (2.8.3) is 
-■'■"•- f f* (2.8.4) 



JO Xi Jo 



/= 



*0 ,2 



2eh xl — z — x 

m 2 2w 2 



To integrate (2.8.4), let 

Z = 
where 



•v-j 



\Zk\x\ - d m y- + dV 



k 2 = * i - 



Vl + B 



(2.8.5) 



(2.8.6) 
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Fig. 2.8.2 — Acceleration-time curves for cushioning with cubic elasticity. See 
equation (2.8.14). 



and B and d m are as given in Part I: 



B = 



4W 2 hr 



d m = J / j/|(-i + vTTi)- 



(1.5.3) 

(1.5.4) 



408 BELL SYSTEM TECHNICAL JOURNAL 

Then (2.8.4) becomes 



- 1 jf 

C0 C JO 



dZ (2.8.7) 



T2 



V(i - 2 2 )(i - £ 2 z 2 )' 

in which the integral is the elliptic integral of the first kind (see Hancock 
"Elliptic Integrals," John Wiley and Sons, New York, 1917). In (2.8.7), 

o, c = coo (1 + B) 1 ' 4 , (2.8.8) 

where co = Vh/nh is the radian frequency that would obtain if the cushion- 
ing were linear with spring rate k . The motion for the linear case has a 
half period, or pulse duration t = tt/coo . The half-period (t 2 ) of the motion 
with cubic elasticity is twice the time required for x 2 to increase from 
to dm. 
From 2.8.5 

lZl -° - °' (2.8.9) 

[Z] X2=dm = 1. 

Hence, from (2.8.7), the half -period is 

= 1 f / % = 2 JE (2.8.10) 

«cA V(l - Z 2 )(l - ^ 2 Z 2 ) C0 C 

where A: is the complete elliptic integral of the first kind. The duration of 
the acceleration pulse is therefore 2K/w c . We can define a radian frequency 
of the acceleration by 

W2 = I = ™ e = Wl + ^) 1/4 (2.8.11) 

0,2 t 2 2K 2K 

The ratio «o/c32 (i.e., t 2 /t ) is plotted in Fig. 2.8.1 which illustrates how the 
pulse duration decreases as the parameter B increases. Hence, for a given 
cushioning with cubic elasticity, the pulse duration decreases as the height 
of drop increases. This is in contrast to the linear case in which the dura- 
tion is independent of the height of drop. 

To find the acceleration x 2 , we return to (2.8.7) and write it in the form 
of an elliptic function: 

snu c t = Z. (2.8.12) 

Substituting the expression for Z given in (2.8.5) and solving for x 2 , we find 

x 2 = d m cn(o> c t - K). (2.8.13) 

Finally, differentiating (2.8.13) twice with respect to /, we find the accelera- 
tion to be 

x 2 = ttUM&snKuct - K) -l]cn(u> c t - K). (2.8.14) 
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The ratio — x 2 /G g is plotted in Fig. 2.8.2 against a radian coordinate 
(a>oO for several values of B. It may be seen that, as B increases, the maxi- 
mum acceleration increases, the duration of the pulse decreases (see Fig. 
2.8.1) and the acceleration-time curve becomes bell shaped. For reference, 
the sinusoid for the linear case (B = 0) is plotted in the figure. 

Figure 2.8.2 is plotted for perfect rebound. If rebound does not occur, 
the curves continue, mirrored in the time axis, so as to form a periodic 
vibration of period 2t2 . 

2.9 Acceleration-Time Relation for Tangent Elasticity 

In this section the effect of tangent elasticity on the shape of the accelera- 
tion-time relation will be studied. The shape of the load displacement 
curve is given by 

2k db . irx 2 . . 

P = tan — . (1.4.3) 

■k Zdb 

The system considered is again that shown in Fig. 1.2.1. Referring to the 
energy equation (1.2.13): 



M 2 X 2 



+ / Pdx 2 = m 2 gh, (1.2.13) 

Ja 



we substitute the above value of P and perform the indicated integration to 
obtain, for the velocity, 



*- A +8 ^ logcos S- (2 - 9 - x) 



Then, as in Section 2.8, 

' 2 dx% f Z2 dx2 



r 2 dx 2 r 2 dx2 

t = I — = / — , = (2 9 2' 

Jo X2 Jo . / n , , 8£ d\ . irx 2 ' ' 



+ 



2f* + =3 log cos =J 

m 2 ir 2 ° 2db 



and the half-period (t 2 ) of the motion is again twice the time required 
for x 2 to increase from to d m . Hence 



r Qm dx 2 

t 2 = 2 / / -, , 2 (2.9.3) 



<A»+^ h «~a 



where, from Section 1.9, 



[-mi 



d m =^cos- 1 ex P | -",-(7) I- :"^-^ 
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The radian frequency of the acceleration is 



7T 

W2 = — 

T2 



and this is to be compared with the frequency 

v _ /fa 
to V tlh 



O>0 — — 



that would obtain if the cushioning were linear with spring rate fa . The 
ratio wo/w2 (i.e., t 2 /t ) was obtained by numerical integration of (2.9.3) 
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Fig. 2.9.1— Duration of acceleration pulse for cushioning with tangent elasticity. See 

equation (2.9.3). 

and is plotted in Fig. 2.9.1 against the ratio d b /d . The figure shows that 
for db/do < 1, the pulse duration varies almost linearly with d b /d . As the 
bottoming distance becomes larger than that required for linear cushioning 
with spring rate fa, the pulse duration approaches asymptotically the 
duration tt/ojo for the linear case. 

As db/do decreases, the pulse duration becomes shorter, but the maximum 
acceleration increases, in accordance with equation (1.9.4) and Fig. 1.9.1. 
The shapes of the acceleration-time curves for several values of db/do are 
illustrated in Fig. 2.9.2. They are more sharply peaked than the corre- 
sponding curves for cubic elasticity (Fig. 2.8.2) as might be expected from 
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the fact that the load-displacement curve for tangent elasticity rises more 
rapidly than that for cubic elasticity; that is, the bottoming is harder. 
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Fig. 2.9.2 — Acceleration-time curves for cushioning with tangent elasticity. 

The curves of Fig. 2.9.2 were obtained by numerical integration of equa- 
tion (2.9.2), to obtain x 2 as a function of /, following which these values 
were substituted in the equation 

.. 2*0 db . 7T.V 2 _ 

m 2 x 2 H tan 7rr =0 

to obtain .i : 2 . It may be observed that the maximum values of the curves 
are the values dictated by equation (1.9.4). 
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In performing the numerical integrations of equations (2.9.2) and (2.9.3), 
it is found that the integrand becomes infinite when x 2 = d m since at this 
point the velocity is zero. In order to avoid this difficulty, it was assumed 1 
that, for a small distance in the neighborhood of d m , the acceleration is 
constant with magnitude G m g as obtained from equation (1.9.4). The 
procedure is described in further detail in Section 2.12. 

Figure 2.9.2 gives the acceleration-time curve for perfect rebound. If 
rebound does not occur, the acceleration is a periodic vibration, each suc- 
cessive half period having the shape shown, with alternating sign. 

2.10 Acceleration-Time Relation for Abrupt Bottoming 

By abrupt bottoming, we mean bilinear cushioning (Class D) as treated 

in Section 1.12. The load-displacement relation is (see equation (1.4.4) 

and Fig. 1.4.4) 

P = k x 2 0^x 2 ^d. 1 

(1.4.4) 

P = k b x 2 — (h — k )db «i > d.j 

Considering, again, the system illustrated in Fig. 1.2.1, the equation of 
motion of m 2 , before bottoming, is 

m 2 x 2 + k x 2 = ^ x 2 ^ d, (2.10.1) 

with initial conditions 

Nm = 0, [* 2 ]m> = ^/2gh. (2-10.2) 

The solution of (2.10.1) is then 

-^S*sin«o*, 0^^„ (2.10.3) 



x 2 = 



where 




av ~~ A/ m 2 



(2.10.4) 



The time (/,) at which x 2 reaches d s is found from (2.10.3): 



1 . _! md, 1 . -id, . 

t, = — sin —nF=i = — sin j , (2.10.5) 

coo V2g« wo do 



where 



-/ ? 



, . /2W 2 h 

da = 



«0 
l 



See Timoshenko, "Vibration problems in Engineering," D. Van Nostrand Co.. New 
York, Second Edition (1937) page 123. 
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i.e., do is the displacement that would have been reached if the spring rate 
remained constant. 
The velocity of m 2 at time /« is 

fcdi-f, = y/2gh cos u t a = A/2gh(\ - -pj. (2.10.6) 

If s/lgh > cood, , the displacement will exceed d„ and the equation of 
motion becomes 

m 2 x 2 + k b x 2 — (k b — ko)d„ = 0, x 2 ^ d, . 
The solution of (2.10.7), with initial conditions 



(2.10.7) 



(2.10.8) 



is 



x 2 = 



where 



kodo 



+ 



+ H-:) rf - 



(2.10.9) 



x 2 ^ </, 



tan" a = 



o>& = 



a-). 



(2.10.10) 



By differentiating (2.10.3) and (2.10.9) twice with respect to /, the 
accelerations for the two regions are found to be 

*t = - G g sin wo/, ^ * 2 ^ <*., (2.10.11) 



•V'2 = — 
#2 ? </«, 



where 



Go = 



</■ 



2hko 



(2.10.12) 



(2.10.13) 



Typical shapes of the acceleration pulse represented by equations (2.10.11) 
and (2.10.12) are shown in Fig. 2.10.1. The curves are drawn for d,/d — 



414 



BELL SYSTEM TECHNICAL JOURNAL 



0.5 and for several values of h/h . The peak values of the curves are the 
same as given by equation (1.12.3). The curve marked k b /k = 1 is the 
sinusoid of the linear case with duration 



(2.10.14) 




Fig. 2.10.1 — Acceleration-time curves for cushioning with bi-linear elasticity. 
d./do = 0.5. See equations (2.10.11) and (2.10.12). 

As before, if the package does not rebound, the acceleration shown is mir- 
rored in the time axis after each half cycle, to form a vibration of period 2r 2 . 

It is useful to know the duration of the complete pulse (aa' in Fig. 2.10.1) 
and also the duration of bottoming (bb' in Fig. 2.10.1). Calling the former 
to and the latter t b , we have, from equations (2.10.11) and (2.10.12) 
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Fig. 2.10.2 — Pulse durations for cushioning with bi-linear elasticity. See equations 

(2.10.15) and (2.10.16). 



(2.10.15) 




(2.10.16) 



These two equations are plotted in Fig. 2.10.2 for several values of k b /k u . 

2.11 Acceleration-Time Relation for Hyperbolic Tangent 

Elasticity 
The relation between acceleration and time for hyperbolic tangent 
elasticity is found by the same procedure that was used for tangent elasticity 
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in Section 2.9. The system considered is that shown in Fig. 1.2.1 and the 
load displacement curve of the cushioning is given by 

P = P tanh -=- . 

-TO 

Substituting the above expression for P in the energy equation (1.2.13), we 
find the velocity to be 



x 2 = Jlgk - ^l log cosh ^? . (2.11.1) 

Then, as before, 

t = f£ (2.11.2) 

Jo %% 

and the half period (t 2 ) of the motion is twice the time required for x 2 to 
increase from to d m , or 

, = 2^, (2.11.3) 

Jo #2 

( 

where, from Section 1.13, 

*-"-ra-*-*-»CSr)- (U3 ' 3) 

The radian frequency of the acceleration is defined as 

«2 = — 

T2 

and this is to be compared with the frequency 

uo = - = A/ — 
to V m 2 

that would obtain if the cushioning had a constant spring rate equal to the 
initial spring rate (£ ) of the hyperbolic tangent cushioning. The ratio 
wo/o>2 (or t 2 /t ) is plotted, in Fig. 2.11.1, against the dimensionless param- 
eter Po/W 2 G Q (see Section 1.13). It may be observed that the pulse 
duration becomes very long when P /W 2 G is small, i.e., when the horizontal 
portion of the load displacement curve (Fig. 1.4.5) comes into play. The 
influence on the shape of the acceleration-time curve is illustrated in Fig. 
2.11.2. The curve marked P /W 2 G — » « is the sinusoid for the linear 
case. For small values of Po/W&o the curve approaches a square wave. 
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Fig. 2.11.1 — Duration of acceleration pulse for cushioning with hyperbolic 
tangent elasticity. 



2.12 Numerical Procedure for Finding Acceleration-Time Relation 
for Class F Cushioning 

When the load-displacement curve does not resemble one of Classes A to 
E, the acceleration-time relation may be found by numerical integration. 
Combining the energy equation, 



nil ±2 



) + J P dx 2 = niigh, 



(1.2.13) 
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with the equation relating time and velocity, 

JO A'2 

we find 



r x -dxi 
" Jo fe' 



t = 



2g Jo 



dxt 



W*h 



-f 

Jo 



Pdx 2 



(2.12.1) 



(2.12.2) 



As an example, consider the problem of a 15 -pound article supported on 
cushioning with the load-displacement curve shown in Fig. 2.12.1. The 
package is to be dropped from a height of 3 feet. The computations are 
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Fig. 2.11.2— Acceleration-time curves for cushioning with hyperbolic tangent elasticity. 

given in detail in Tables III and IV. The headings of Columns (1) to (8) 
of Table III are the same as in Table II, Section 1.15. A n is the integral 
under the radical of equation (2.12.2). Column (10) of Table III is the 
integrand of Equation (2.12.2), i.e., it is proportional to the reciprocal 
of the velocity expressed as a function of displacement. The function 
is plotted in Fig. 2.12.2 and its integration is performed in Table IV. In 
columns (11), (12) and (13), intervals of x 2 are chosen to suit the shape of 
the curve. The values for column (14) are taken from column (10). Col- 
umns (15) and (16) perform the same operations on the integrand 
(W z h - A„)~ h that are performed in Columns (5) and (6) of Table III on 
the integrand P. 
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Fig. 2.12.1 — A load displacement curve for Class F cushioning. 
TABLE III 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




A(*i)„ 


(*«)» 


^CF.+*.-d 





•.-fe 


*-fc 


II V/ - A n 


l 




VWth - An 


n 





0.0 














540 


0.0431 


i 


.20 


0.2 


105 


10.5 


10.5 


0.7 


7.0 


529.5 


0.0435 


2 


.20 


0,4 


155 


26.0 


36.5 


2.4 


10.3 


503.5 


0.0446 


3 


.20 


0.6 


192 


34.7 


71.2 


4.8 


12.8 


468.8 


0.0462 


4 


.20 


0.8 


217 


40.9 


112.1 


7.5 


14.5 


427.9 


0.0483 


5 


.20 


1 


237 


45.4 


157.5 


10.5 


15.8 


382.5 


0.0511 


6 


.20 


1.2 


257 


49.4 


206.9 


13.9 


17.1 


333.1 


0.0547 


7 


.20 


1 4 


277 


52.9 


259.8 


17.3 


18.5 


280.2 


0.0597 


8 


.20 


1.6 


305 


58.2 


318.0 


21.2 


20.3 


222.0 


0.0671 


9 


.20 


1.8 


342 


64.7 


382.7 


25.5 


22.8 


157.3 


0.0798 


10 


.20 


2.0 


392 


73.4 


456.1 


30.4 


26.1 


83.9 


0.109 


11 


.05 


2.05 


405 


19.9 


476.0 


31.8 


27.0 


64.0 


0.125 


12 


.05 


2.10 


422 


20.7 


496.7 


33.2 


28.1 


43.3 


0.152 


13 


.05 


2.15 


440 


21.6 


518.3 


34.6 


29.4 


21.7 


0.215 


14 


.01 


2.16 


445 


4.42 


522.7 


34.8 


29.7 


17.3 


0.240 


15 


.01 


2.17 


450 


4.48 


527.2 


35.2 


30.0 


12.8 


0.279 


16 


.01 


2.18 


455 


4.52 


531.7 


35.5 


30.3 


8.3 


0.347 


17 


.01 


2.19 


457 


4.56 


536.3 


35.8 


30.5 


3.7 


0.521 


18 


.01 


2.20 


462 


4.60 


540.9 


36.1 


30.8 





oo 
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TABLE IV 



(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 




A(»)» 


(*s) n 


/» = 
1 
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Fig. 2.12.2— Plot of Column (3) vs. Column (10) of Table III. 

A difficulty arises because the integrand (W 2 h — ^4„) -i becomes infinite 
for the maximum displacement (see Column (14)). This is avoided by 
assuming that the acceleration is constant in the last interval 2 and has the 

2 Timoshenko, "Vibration Problems in Engineering," D. Van Nostrand Co., New York, 
Second Edition (1937) page 123. 
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value given in Column (8), Table III, for the maximum height of drop. 
Then, 

A(* 2 ) n = hG ng At> (2.12.3) 

or 



A* = 



V 



2A(*2) n 
G m g 



(2.12.4) 
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Fig. 2.12.3 — Acceleration-time curve (for the cushioning shown in Fig. 2.12.1) obtained 

by numerical integration. 

In the present instance, 

(Ax 2 ) n = 0.01 inches 

G m g = 30.8 X 386 = 11900 in/sec. 2 

Hence, from (2.12.4), A/ = 0.0013 sec. and the last entry in Column (17) 
is obtained by adding this value of A/ to the preceding entry. 

The final curve of acceleration vs. time is obtained by plotting the entries 
of Column (17) against the entries of Column (8), Table III, for correspond- 
ing values of x% . The result is shown in Fig. 2.12.3. 
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PART III 

AMPLIFICATION FACTOR 

3.1 Introduction 

If the maximum acceleration, of the packaged article as a whole, is 
reached very slowly, the severity of the disturbance experienced by a 
structural element of the packaged article is very nearly proportional to the 
maximum acceleration. Roughly speaking, "very slowly" means that the 
time, during which the acceleration undergoes a major change in magnitude, 
is long in comparison with the natural period of vibration of the element 
under consideration. When this is so, no transient vibration is excited in 
the element. The displacement response of an element under very slowly 
varying conditions is called the "static response". Under more rapidly 
varying conditions the dynamic response to the same maximum acceleration 
may be greater or less than the static response. The ratio (A) of the maxi- 
mum dynamic response to the static response is called the amplification 
factor. In general, for a given acceleration disturbance, very low-frequency 
elements have amplification factors less than unity, while the amplification 
factors are greater than unity for elements whose natural frequencies are 
near or above the disturbing frequencies. The numerical value of the 
amplification factor depends not only on the manner in which the disturbing 
acceleration varies with time, but also on the "reference acceleration", i.e., 
the value of acceleration for which the static response is calculated. Usually 
the reference acceleration chosen for calculating the static response is the 
maximum value (G m ) of the disturbing acceleration. However, when 
special circumstances are being investigated, such as the effect of damping 
or abrupt bottoming, the reference acceleration is taken to be Go , which is 
the acceleration that would be reached if the damping or bottoming were 
absent. In such cases the amplification factor includes both the effect of 
rate of change of acceleration and the effect of the special conditions. 

When the reference acceleration is G m the amplification factor will be 
denoted by A m and when the reference acceleration is G the amplification 
factor will be denoted by A . The symbol G e will be used to designate the 
slowly applied acceleration that would produce the same maximum dis- 
placement as the transient acceleration, i.e., G e = A m G m or G e = A(P • 
The symbol G 3 will be used to denote the safe value of G e , for an element 
of the packaged article, as determined by a strength test or by calculation. 
In specifying G„ some judgement is required to take into account the effects 
of plastic deformation in comparing tests made on greatly different time 
scales. Good judgement is also necessary in deciding whether or not the 
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assumptions listed in Section 0.2 are valid in each application. The general 
procedure for using amplification factors is as follows. We first find the 
value of the reference acceleration (in units of number of times gravity) from 
Part I. From Part II we find the properties of the acceleration-time rela- 
tion which give us the information required for entering one of the curves 
of Part III and finding the amplification factor. Then, the product of the 
reference acceleration and the amplification factor (A m G m or A^Go) is a 
number (G e ) by which the weight of the structure is to be multiplied when 
calculating its deflection or stress by the usual static methods of elementary 
strength of materials. Alternatively, G, must be found not to exceed G s . 




i- -J: 



(a) 



(b) 



Fig. 3.2.1— Idealized system used in calculating amplification factors for linear undamped 
cushioning with perfect rebound, (a) initial position, (b) first contact with floor. 

In the following sections the amplification factors for typical transient 
accelerations encountered in package drop tests are calculated. The ampli- 
fication factor curves that are plotted are entirely analogous to the familiar 
"resonance curves" for steady sinusoidal vibration, except that in this 
case the disturbing forces are transients of various shapes. It will be seen 
from the curves that the maximum acceleration, as calculated by the 
methods of Part I or as measured by an accelerometer, is not always a true 
measure of the severity of the disturbance. 



3.2 Amplification Factors for a Half-Sine-Wave Pulse 
Acceleration 

The first case to be treated is the response of an element of the packaged 
item to the transient acceleration that would occur in a package with linear 
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undamped cushioning and perfect rebound. Figure 3.2.1 illustrates the 
idealized system, and it may be noted that the mass m 3 is omitted, as is 
required for perfect rebound (Section 2.3). At first we shall consider that 
the mass mi is undamped and later we shall consider the effect of damping 
in this element. 

The mass mi is taken to be small in comparison with m 2 , so that the motion 
of the latter is the same as we found it to be in Section 2.2 where m x was not 
considered. Hence the acceleration of m 2 is a half -sine wave pulse: 

x 2 = —co 2 \/2ghswwst, (0 ^ t ^ 7r/w 2 ). (3.2.1) 

The equation of motion of mi is 

Witfi + *i(*i - »%) = 0. (3-2.2) 

Let x be the relative displacement of mi with respect to m 2 , i.e., 

x = Xi — %$ . (3.2.3) 

x is proportional to the force in the spring (fax) and to the acceleration 
of i»i and hence is proportional to the deflection, strain and stress in the 
element which the system Wi , k x represents. 
Substituting (3.2.3) in (3.2.2), we find: 

mix -f fax = -miXi . (3.2.4) 

This equation holds for the duration tt/co,. of the pulse x 2 . The initial 
conditions for x are 

[x] t=0 = [*] M = (3.2.5) 

so that the solution of (3.2.4) is 

x = V^M [<* sin aM - sin «,*! (o ^ t^ ?). (3.2.6) 

It may be seen that x is composed of a forced displacement at the accelera- 
tion frequency co 2 , on which is superposed a free vibration at the natural 
frequency, m , of the element. The maximum value of the relative displace- 
ment is 



Smax = —. ' v Sin 



Wi 



— ^ Olll 

(* - A s 

\0)2 / W2 



max 

+ 1 



( 0<t< ^)> (3.2.7) 



in which n is a positive integer chosen so as to make the sine term as large 
as possible while the argument remains less than ir. 

(3.2.7) gives the maximum dynamic response of the element Wi during 
the interval of impact. To find the amplification factor we must compare 
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x max with the "static response" i.e. with the value (x 8t ) that x would have 
if the acceleration x 2 reached the same maximum value («2\/2gA) in a very 
long time. The resulting value may be found from (3.2.4) by omitting 
the transient term mix. Then 

T 

Xst = C02 \/2gh — 
or 

OM / 

x lt = -2\/2^h. (3.2.8) 

Oil 

The amplification factor for the interval ^ / ^ 7r/o) 2 is then 

Am m *= . _J*_ sin J!^_ , (O^/^ »/<*)• (3.2.9) 

*.< Wl _ 1 ^4-1 

0)2 0)2 

It should be observed that A m depends only on the frequency ratio a>i/o>2 . 
That is, since o>i/a> 2 = ti/ti , the amplification factor depends only on the 
ratio of the pulse duration to the half period of vibration of the element. 

Thus far we have studied only the motion in the interval ^ / ^ tt/o)2 . 
We must not, however, overlook the possibility of larger displacements of Wi 
with respect to m« occurring after rebound. In fact, examination of (3.2.6) 
reveals that x has no maximum in the interval ^ t ^ 7r/o>2 when on < 0)2 . 
It is very likely, then, that larger values will occur at later times. 

After rebound, Wi executes free vibrations with respect to m* . We have 
to compare the magnitude of .r max , in the interval < / ^ tt/o)2 , with the 
amplitude of the free vibration. Calling the relative displacement during 
free vibration x' and measuring a time coordinate /' from the instant the 
package leaves the floor, we have 

mix' + fax' = 0, (3.2.10) 



with initial conditions 



M«'=0 = M<=r/u, , 



(3.2.11) 



\X IfmO = [X\t=r/o, 2 • 

The solution of (3.2.10) with initial conditions (3.2.11) is 



*' = 



y 4 K 1 + C0 O) • ( ^w\ 

-i j\ ^ '- smlcoi/' + r— J, 

o)i(o>;> — (oi) \ /u>2 / 



(3.2.12) 



/ = * 

0)2 
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Then 
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A m = 



x at 



„OJl 7TC0, 

2-cos — 

0)2 iiC02 



\ -> ««/' 



(3.2.13) 



We find, on comparing (3.2.13) with (3.2.9) that for wi < o> 2 equation 
(3.2.13) gives the larger value of A m , while for on > o> 2 equation (3.2.9) 




Fig. 3.2.2— Amplification factors for linear undamped cushioning with perfect rebound. 
See Fig. 3.2.1 and equations (3.2.9) and (3.2.13). 

gives the larger value of A m . That is, when the duration of impact is 
shorter than the half-period of vibration of the element, the maximum 
displacement (and stress) in the element occurs after the impact is over. 

The curve marked ft = in Fig. 3.2.2 is a plot of the largest value of A m 
from (3.2.9) and (3.2.13) with the frequency ratiowi/co 2 as abscissa. (3.2.13) 
was used for wi/w 2 ^ 1 and (3.2.9) for wi/«a ^ 1. The maximum value 
of A m is 1.76 and occurs atco 1 /co 2 = 1.6. Hence, at this frequency ratio, the 
deformation of the element is 1.76 times as great as would be expected from 
a calculation using the maximum value of acceleration alone as in Part I. 
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On the other hand, for frequency ratios wi/w 2 < 0.5 the severity of the 
shock can be very much less than might be expected from the calculations 
of Part I. For very small values of w 1 /oj 2 the amplification factor may be seen 
from (3.2.13) to be equal to 2a>i/a> 2 . For large values of wi/coo (stiff elements) 
Fig. 3.2.2 shows that the amplification factor is very nearly unity and the 
methods of Part I can be used without additional calculation. 

When damping of the element of the packaged article is considered, the 
amplification factors are less than without damping. The applicable 
equations of motion during and after impact are obtained by inserting 
velocity damping terms in (3.2.4) and (3.2.10): 

mix + ci± + kix = -mix 2 , ^ / ^ — (3.2.14) 

(02 

mix' + cix' + kix' =0, t ^ — . (3.2.15) 

0)2 

If we express the damping of the element mi as the fraction of critical 
damping 

(as in Section 2.5) equations (3.2.14) and (3.2.15) become 

x + 2/3i wi* + <£x = -ft , ^ t ^ — , (3.2.17) 

0)2 

*' + 2/3 1 co l .v' + u\x' =0, t% — . (3.2.18) 

a>2 

The amplification factors for equations (3.2.17) and (3.2.18), with boundary 
conditions (3.2.5) and (3.2.11), respectively, were obtained on the Westing- 
house Mechanical Transients Analyzer 3 iorfii = 0.005, 0.01, 0.05, 0.10, 0.30, 
0.50 and 1.00. The curves are shown in Fig. 3.2.2. 

3.3 Application of Half-Sine-Wave Amplification Factors 

As an example of the use of the amplification factor curves of Fig. 3.2.2, 
let us consider the following problem: 

3 Arrangements for performing these calculations were made through the courtesy of 
Mr. A. C. Monteith, Manager of Industry Engineering, and Mr. C. F. Wagner, Manager 
of Central Station Engineering, Westinghouse Electric and Manufacturing Co. Dr. G. D. 
McCann, Transmission Engineer, was in immediate charge of the project. For a descrip- 
tion of the analyzer see "A New Device for the Solution of Transient- Vibration Problems 
by the Method of Electrical-Mechanical Analogy" by H. E. Criner, G. D. McCann and 
C. E. Warren, Journal of Applied Mechanics, Vol. 12, No. 3 (1945) pp. A-135 to A-141. 
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It is required to judge the suitability of a proposed package for a large 
vacuum tube weighing 10 pounds. Strength tests have been made on the 
tube in a shock testing machine which produces a half-sine-wave acceleration 
pulse of 25 milliseconds duration. The weakest element of the tube is 
found to be the cathode structure, for which the safe maximum acceleration 
in the drop testing machine is 200g. The cathode structure has a natural 
vibration frequency of 120 cycles per second and has 1% of critical damping. 
The proposed package has essentially linear, undamped cushioning with a 
spring rate of 3300 pounds per inch and an available displacement of £ inch. 
The outer container weighs much less than the tube so that the package may 
be expected to rebound. Is the cushioning suitable for protecting the 
cathode in a drop of 5 feet? 

First find the maximum G that the tube will experience in a 5 ft. drop of 
the package (equation 1.3.3): 

r A /2*fe _ . /2 X 60 X 3300 __ lgg 

Gm ~ y w* V io 

The accompanying maximum displacement is, from equation (1.3.4), 

2h 2 X 60 n , . 

d m = — = = 0.6 m. 

m G m 199 

The available displacement (f inches) is therefore sufficient and the maxi- 
mum acceleration (199g) is slightly less than the safe maximum acceleration 
(200g) found with the shock testing machine. However, before the cush- 
ioning is approved it is necessary to investigate the frequency effects. The 
duration of acceleration in both the shock machine and in the package 
must be considered. 

The amplification factor for the element tested in the shock machine is 
found as follows. First find the frequency corresponding to the 25 milli- 
second pulse: 

J2 2 X .025 F 

The ratio of the element frequency to the shock machine frequency is 

fj = wj = 120 = 6 

/ 2 C02 20 

Entering Fig. 3.2.2 witho>i/a) 2 = 6, we read, from the curve fa = 0.01, A m = 
1.14. The 200g test in the shock machine is, therefore, equivalent to a 
slowly applied acceleration of G„ = 200 X 1.14 = 22Sg. 
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To find the corresponding quantity for the package drop, first find the 
cushion frequency: 

/ 1 "j/5 1 A /3300 X 386 M 
The ratio of the element frequency to the package frequency is therefore 

A = «i m 120 ^ 2 x 

ji o)2 57 

Entering Fig. 3.2.2 with o>i/co2 = 2.1 we read, from the curve ft = 0.01, 
A m = 1.59. The 199g acceleration pulse in the package drop is therefore 
equivalent to a slowly applied acceleration of G e = 199 X 1.59 = 316g. 
This is almost 40% in excess of the value (228g) found to be safe from the 
shock machine data. The cushioning is therefore judged to be inadequate. 
Tne procedure for finding the correct spring rate for the cushioning is as 
follows. It is known that we must have 



Therefore, take 



Now 



Therefore 



Also 



G e ^ G s . 



G e = A n G m = 228. 



W 2 t = V!' 



A m u 2 = 409 rad/sec. 



«i = 2x X 120 = 754 rad/sec. 

Then, with successive trial values of w 2 , we calculate coi/wo , enter Fig. 3.2.2, 
read the corresponding value of A m from curve /3i = 0.01 and test to see if 
the product A m u)2 = 409. The combination which satisfies the test is found 
to be 

w 2 = 280 rad/sec. 

ioi/io-t = 2.69 

A m = 1.47. 



Then 
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h » £ m = gggy< 10 = 2030 lbs./in. 



-^ 



'■■- - ^i? ■ 155 

*. = £ - .77 in. 

Hence the spring rate of the cushioning should be reduced from 3300 lbs./in. 
to 2030 lbs./in. and the available space should be increased to accomodate 
the 0.77 inch maximum displacement before bottoming. 

3.4 Special Treatment of Strong, Low Frequency Elements 
The product of the amplification factor (A m ) and the maximum accelera- 
tion (G m ) must be equal to or less than the maximum allowable slowly 
applied acceleration (G„): 

G e = A m G m ^ G, . (3.4.1) 

For frequency ratios 

^<i (3.4.2) 

Figure 3.2.2 shows that, approximately, 

A m = 2^ (3.4.3) 

for Pi ^ 0.10. When this is so, we may combine (3.4.1) and (3.4.3) to 
obtain the criterion 



2^G m ^G.. . (3.4.4) 

Oil 



Now, 



Hence the criterion (3.4.4) may be written as 

2*>i J- ^ G s (3.4.6) 



or 



where h is in inches and /i is in cycles per second. 



A 3 ^7= (3A7) 
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It may be observed that (3.4.7) is independent of the properties of the 
cushioning. Hence, as long as (3.4.2) is satisfied, any cushioning at all 
may be used for an element that satisfies (3.4.7) regardless of the magnitude 
of the maximum acceleration G m . In particular, rigid mounting is suitable 
for such an element. The only precaution to be observed is that the maxi- 
mum acceleration and duration must not be unfavorable for other elements 
of the packaged article. 

Example: A 9-pound vacuum tube has an anode structure for which the 
safe maximum acceleration is 200g as determined in a centrifuge test. 
The natural vibration frequency of the anode is 35 cycles per second and the 
damping is 1% of critical. What cushioning around the tube is required to 
protect the anode from damage in a package drop of 3 feet? 

Calculate 



1.1G„ 
Vh 



1.1 X 200 

~vw = 36 - 7 cps - 



This is greater than f\ = 35 c.p.s and hence any cushioning is safe for the 
anode. The results of calculations for cushioning with spring rates of 50, 
500, 5000, 5 X 10 and 5 X 10 pounds per inch are given in the following 
table: 



*i(lbs./in.) 


G m 


- 


A m 


A m G m 


50 


20 


4.74 


1.12 


22 


500 


63 


1.47 


1.65 


106 


5 X 10 3 


200 


.474 


0.9 


180 


5 X 10 5 


2,000 


.0474 


0.09 


180 


5 X 10 7 


20,000 


.0047 


0.009 


180 



In each case the product of A m G m is less than the allowable 200 and, as long 
as the combination of G m and the amplification factors for other elements 
does not exceed the allowable A m G m for those elements, the cushioning is 
suitable. The precaution to observe is that higher-frequency elements 
shall not have amplification factors such that A„,G m may be excessive for 
them. 



3.5 Amplification Factors for Damped Sinusoidal Acceleration 

If the outer container of the package is heavy enough (see Section 2.3) 
there will be no rebound and the packaged item will vibrate in the cushion- 
ing after impact. For linear cushioning with velocity damping, the accelera- 
tion produced by the vibration will be a damped sinusoid (equation (2.5.5) 
and Fig. 2.5.2). The system to be considered is shown in Fig. 3.5.1. 
To determine the effect of the damped vibration of ra 2 on the mass »j, , 
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(a) 

Fig. 3.5.1— Idealized system for linear damped cushioning with no rebound, 
(a) initial position, (b) first contact with floor. 
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Fig. 3.5.2 — Amplification factors for linear damped cushioning with no rebound. 
/3 2 = 0.005. See equations (3.5.1) and (3.5.2). 

we note that the equation of motion and initial conditions are identical 
with (3.2.17) except that for the acceleration x 2 we use the damped sinusoid, 
equation (2.5.5), instead of the half-sine pulse (3.2.1). The solution of 
(3.2.17), i.e. the relative displacement (*i - x 2 ) of wi with respect to m 2 , 
is 
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= 4^ le^'U sin (di + 7 - a) + B sin («[/ - 7 - fll 
- e _ ^ 2 'U sin («{/ + 7 - 5)- 5 sin (co 2 * + y + f )] } (3.5.1.) 



where 







. / 


'h 






Oil 










(i)o 


— A/ 








= 


i 

WlVl 

oj* Vl 


m% 




COi 


— 


0? 


CD] 




01 



l/.'l = \/(/3 2 w 2 - /3io>i) 2 + («{ - «s)= 

1/5 = V(02CO2 - ft Wl ) 2 + («{ + O)0 2 

2$ - 1 

tan T = 2* vi - 3 

«1 — 0)2 
tan 6 = 



tan f = 



02 0)2 — 01 0>i 
0)1 ■+■ 0)2 



?2 0)2 — ftWl 

The relative displacement of m% with respect to m% is seen to consist of a 
forced, damped vibration (o) 2 , /3 2 ) on which is superposed the free damped 
oscillations (0)1 , /3i) of m\ . 

The amplification factor 

Ao . *- _ ^^ (3.5.2) 

is plotted in Figs. 3.5.2 to 3.5.7 for six values of jSi and six values of /3 2 . 
These curves were obtained by direct solution of the differential equation 
on the Westinghouse Mechanical Transients Analyzer. The amplifica- 
tion factor in this case includes the effect of damping; i.e., the reference 
acceleration is not the maximum acceleration of mi , but is the maximum 
acceleration that nii would reach if the damping /3 2 were zero. Conse- 
quently, the amplification factors for large values of o)i/o>2 do not approach 

4 See footnote, Section 3.2. Only enough data were obtained with the analyzer to find 
the general shapes of the curves, so that the fine structure is not revealed. Checks on 
the analyzer results were made by computing A§ from equations (3.5.1) and (3.5.2) for 
W1/0J2 = 1, ft = /3 2 ; wi/wj = 0; a>i/a> 2 — » co . 
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Fig. 3.5.3— Amplification factors for linear damped cushioning with no rebound, 
fo = 0.01. See equations (3.5.1) and (3.5.2). 



I 

z 


H 
< 




Fig. 3.5.4— Amplification factors for linear damped cushioning with no rebound. 
Pz = 0.05. See equations (3.5.1) and (3.5.2). 

unity. For example, the curve for ft = 0.005, 0, = 1 (Fig. 3.5.7) approaches 
a value of nearly four as a>i/co 2 -» °° • The factor four is composed of two 
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Fig. 3.5.5 — Amplification factors for linear damped cushioning with no rebound. 
fc = 0.1. See equations (3.5.1) and (3.5.2). 
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Fig. 3.5.6 — Amplification factors for linear damped cushioning with no rebound. 
2 = 0.5. See equations (3.5.1) and (3.5.2). 

factors of two. The first arises from the fact that the maximum value of 
acceleration, for /3 2 = 1, is twice the value that would be reached if /3 2 
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were equal to zero (see Fig. 2.5.3). The second factor (of nearly two) is 
due to the fact that the maximum acceleration is reached at time t = 
when /3 2 = 1 (see Fig. 2.5.2) and the response of an almost undamped sys- 
tem (fii = 0.005) to a suddenly applied and subsequently maintained 
acceleration is double the response to a slowly applied acceleration (see 
curve (a) Fig. 3.8.1). For /3i > and /3 2 < 1 the amplification factor is 
less than four, as «i/o> 2 — > °° , in accordance with the curves plotted in Fig. 
3.5.8. 

Example: A 1.5-pound vacuum tube is to be packed in a container whose 
estimated weight will be at least 50 pounds. The cathode structure of the 
tube has a natural frequency of 25 c.p.s. with damping 0.5% of critical 
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Fig. 3.5.7— Amplification factors for linear damped cushioning with no rebound. 
/3 2 = 1.0. See equations (3.5.1) and (3.5.2). 

and its safe acceleration, as determined in a centrifuge, is 90g. What 
spring rate of cushioning is suitable for protecting the cathode in a drop of 
five feet? It is specified that the cushioning shall have damping 50% of 
critical. 

Assuming linear cushioning, the spring rate that would be prescribed, by 
considering maximum acceleration alone, is 

g 2 g > = 1.5X(90) 2 = 1011bs/ . n 
' 2h 2 X 60 

Considering damping, Fig. 2.5.3 shows that 50% of critical damping does 
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not change G m . To find the amplification factor we must first decide if the 
package will rebound. With 50% of critical damping, the maximum ac- 
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Fig. 3.5.8 — Limiting values of amplification factors for linear damped cushioning with no 
rebound. coi/w 2 — > oo. See equations (3.5.1) and (3.5.2). 

celeration on the first upstroke is 0.164 G m (see Section 2.6 and Fig. 2.6.1). 
Then, 0.164 X 90 X 1.5 = 22 lbs. which is less than the estimated weight 
of the outer container. The package will not rebound and Fig, 3.5.6 
should be used for the amplification factor. 
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The frequency of vibration of the tube in its cushion will be 

Hence Wl /co 2 = 2tt X 25/159 = 0.99 and, from Fig. 3.5.6, A = 1.4. Hence 
G e = 90 X 1.4 = 126, which is greater than the allowable G, = 90, so 
that the 101 lb ./in. cushion is unsatisfactory. 
To obtain satisfactory cushioning, set 

AoG = 90, 



that is 



90 
A o>2 = —/=i = 159 rad/sec. 



Noting thatwi = 2tt X 25 = 157, we find from Fig. 3.5.6 that there are two 
values of w 2 (90 and 600 rad/sec.) that satisfy the criterion A^ 2 = 159 
rad/sec. The first gives 

co 2 = 90 rad/sec. 

^o= 1-8 

k 2 = 31.5 lbs./in. 

Go = 50 

G e = 90 

d m = 2.4 in. 



The second gives 



c>2 = 600 rad/sec. 

A = 0.27 

ko = 1400 lbs./in. 

Go = 335 

G e = 90 

d m = 0.36 in. 



The second solution requires less space for cushioning than the first but 
should be used only if the remainder of the tube can endure the high ac- 
celeration of 335g. Otherwise the 50g package should be used . 
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3.6 Amplification Factors for the Pulse Acceleration of 
Cubic Cushioning 

In a rebounding package with undamped Class B cushioning, the pack- 
aged article (»i 2 ) will undergo a pulse acceleration of duration tt/wz as given 
by equation (2.8.11). The shape of the pulse is illustrated in Fig. 2.8.2 
and its functional form is 

,, = *J^ [tf rf (25*' - *) - l] «< (^' - K) . (2.8.14) 

To determine the influence of the shape and duration of this pulse on the 
amplification factor, we proceed as before by substituting (2.8.14) in the 



S 




///////////////////// 



7777? 



77777777777777777 



(a) (b) 

Fig. 3.6.1— Idealized system used in calculating amplification factors for non-linear, 
undamped cushioning with perfect rebound. 

differential equation governing the relative displacement (x = Xi — x 2 ) 
between m\ and w 2 (see Fig. 3.6.1): 

x +ta\x = -xi. (3.6.1) 

With boundary conditions .t(0) = x(0) = 0, the solution of (3.6.1) may be 
written as 

1 f' 
x as — I x 2 (\) sin wi (X — /) 
o>i Jo 

and the maximum value of * may be expressed by 

1 r' m 

x m ax *= — / :v 2 (X) sin o>i(X — t m ) dX, 
where l m is the time at which the largest value of x occurs. 



dX 



(3.6.2) 



(3.6.3) 
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The amplification factor, in this case, will be taken as the ratio of x m „ 
to the relative displacement (x st ) resulting from a slow application of the 
maximum value of .f 2 . From (3.6.1), 



X sl = 



X-2 



2 






where G m is given by equation (1.5.6). Then 



A m = 



x al 



G m g 



(3.6.4) 



(3.6.5) 





Fig. 3.6.2 — Amplification factors for undamped cushioning with cubic elasticity. Perfect 
rebound. See equation (3.6.6). 



or 



A m = 



G m g 



rim 

\ x 2 {\) sin «i(X - / m ) d\. (3.6.6) 



A m was evaluated, mostly by graphical methods, for four values of B 
(0, 2, 20 and °°) and the results are plotted in Fig. 3.6.2. Observing that 
B = corresponds to linear cushioning, it may be noted that cubic non- 
linearity in the cushioning does not change the amplification factor by 
more than 35% even in the most extreme case (B — > °o). The severity 
of the shock, however, may be much greater for the cubic cushioning than 
for linear cushioning with a spring rate equal to the initial spring rate (£ ) 
of the cubic cushioning. This is because A m is multiplied by G m to obtain 
G e and, for large values of B, G m may be much larger than the maximum 
acceleration for the linear case. In other words, in comparing Class B 
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with Class A cushioning the difference in maximum acceleration, rather than 
the difference in amplification factors, is usually more important. 

Example: Consider the example given in Section 1.6 and let it be required 
to determine the effect of pulse duration on a cathode structure with a 200 
c.p.s. natural frequency of vibration. In Section 1.6 we found that 

B = 5.4 ' k = 255 

Go = 28.6 r = 108. 

G m = 55 

With B = 5.4, enter Fig. 2.8.2 and find 



^ = 0.88. 



w 2 
Now 



Vg-l/^- 4 "- v~ 



co 



Hence 



wa = jj-jjg = 75 rad./sec. 



Then, with wi/o> 2 = 2tt X 200/75 = 16.7, enter Fig. 3.6.2 and find A m = 
approximately 1.0. Hence G e is about the same as G m and the conclusions 
reached for this problem in Section 1.6 are not altered. 

3.7 Amplification Factors for Abrupt Bottoming 

The amplification factors for bilinear elasticity have not been computed 
in complete detail. They can be obtained approximately by using the dura- 
tion curves (Fig. 2.10.2) and the amplification curves for the linear case 
(Figs. 3.2.2 and 3.5.2 to 3.5.7). It is useful, however, to calculate the am- 
plification factors for extremely abrupt bottoming (kb —> °°) to obtain a 
general understanding of the accompanying phenomena. 

The system to be considered is illustrated in Fig. 3.7.1. It is assumed 
that the impact between m 2 and the base (occurring at / = /, , Xj = d,) 
has a coefficient of restitution of unity. Hence m- 2 will strike the base with 
velocity 



r*J«. = iA*(i-g) 



(see equation (2.10.8)) and leave it at a velocity of the same magnitude but 
opposite sign. Perfect rebound of the whole package is also assumed. 
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The acceleration pulse will then look like the curve marked k b /k -* =° 
in Fig. 2.10.1. 

There will be three regions in which to consider the relative displacement x: 

Region 1 < t < t, 

Region 2 /. < t < 2t s 

Region 3 I > 2t. 

The relative displacement (x = x x - x 2 ) of m x with respect to m 2 will have 




Fig. 3.7.1 — Idealized system representing abrupt bottoming. 

the same functional form for Region 1 as in the linear case (see equation 
(3.2.6)), and the amplification factor is, by analogy with (3.2.9), 



A = 



Wo 



sin 



2mr 



55 - 1 55 + 1 

O)0 k-'u 



< / < t, , 



(3.7.1) 



where 



wo = ko/niz . 
For Region 2, we use the differential equation 

x + o>i .t = o>o V 2gA sin aj (/ — 2/,) 



(3.7.2) 



and, as initial conditions at / = /, , we use the terminal conditions for 
Region 1 with the sign of [x 2 ]t=t, reversed. The amplification factor for 
this region is found to be (by the same method as in Section 3.2): 
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WT 



A ° = \ZT1 ^ A * ~*~ B ~ sin ^ Wllm ~*~ v ~ Wl '*^ 



1 . / . _ x d\ (3.7. 

7 \>> sin I ojo /,„ — coo / s — sin -=- J , 

i _ [55 j V V 



/. < t ' < it, 

where 

Wo 



h)0 



jd . "i . /«i • -i d s \ 

7z=t /i = — 7 \ 2 sin I — sin — J 

V2gh i _ ( ^ ] v°° "°/ 



m 



J!!l R w, fo 5i , /i <** fa • -J <*A1 



. -i^i 

»7 = tan -g 

and /„, is the root of 

2 \ZT~h ** A2 ~*~ B2 cos ' W| /m ~*~ v ~ Wl '^ 



"2 cos ( coo t m — wo t a — sin * T-" J = 



" fe )' 



. -i rf.\ 
5m rfoJ = 



that yields the largest value of A in equation (3.7.3). Region 3 is gov- 
erned by 

* + uix = (3.7.4) 

and the initial conditions are the terminal conditions of Region 2. By 
the same method as was used in Section 3.2, we find 



2"i 

Wo 

-In = 



\Wo/ 



r(-y i / 1 _|_ cos (- sin -.Ai, 

! LW r "5 \<»o oo/J (3.7.5 



t> 2t,. 



The largest value of /1 from equations (3.7.1), (3.7.3) and (3.7.5) is 
plotted against coi/w in Fig. 3.7.2 for several values of d,/do . 
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Notice that the amplification factor is A rather than A m . That is, 
the reference acceleration is G rather than G m . This is necessary because 
G m is infinite in the present instance. Hence Fig. 3.7.2 cannot be com- 
pared directly with Figs. 3.2.2 and 3.6.2. However, it is interesting to 
observe that, for «i/w 2 < 0.5, (low frequency elements) abrupt bottoming 
has no harmful effect. For high-frequency elements, the severity of bottom- 
ing is very great even when very nearly all of the required space (d ) is 
available. For example, if 90% of the required space is available (d a /d = 
0.9) and the frequency of the element is ten times the package frequency, 




Fig. 3.7.2— Amplification factors for abrupt bottoming. See equations (3.7.1), (3.7.3) 

and (3.7.5). 

the severity of the shock is almost ten times as great as it would be if the 
additional 10% of space were available. 

3.8 General Influence of Shape of Acceleration-Time Curve 
on Amplification Factor 

When amplification factor curves are not available for a special shape of 
acceleration-time curve, an approximate value of A m may be obtained by 
interpolation between or extrapolation from the curves of the preceding 
sections. The shape of the acceleration-time curve and its duration (t 2 ) 
or frequency (w 2 ) should be found, first, by the methods described in Part II. 
The shape found should then be compared with the standard shapes shown 
in Part II, for which amplification factors are given in Part III. 

The amplification factor found in this way will generally be within 25% 
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of the true value because amplification curves for pulse accelerations do not 
differ greatly even for very different acceleration-time curves as long as the 




Fig. 3.8.1 — Dependence of amplification factor on shape of symmetrical acceleration pulse. 




Fig. 3.8.2 — Effect of asymmetry of an acceleration pulse on amplification factor. 

amplitudes and frequencies are adjusted to the same scales. This is illus- 
trated in Fig. 3.8.1 where the amplification factor curves are drawn for 
square wave, half-sine wave, triangular and cubic pulses. 
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Amplification factors for small values of w 1 /w 2 may be calculated very 
accurately if it is observed that the initial slope of the amplification factor 
curve for a pulse acceleration is proportional to the area under the accelera- 
tion-time curve. For example, noting that the initial slope of the amplifica- 
tion factor curve for the half-sine wave pulse is 2, we assign the value 2 to 
the area under the half-sine wave. On the same scale, the area under a 
square wave pulse is -k and under a triangular pulse is t/2. Accordingly, 
the initial slopes of the amplification factor curves for the latter two pulses 
are ir and x/2 respectively. 

As an additional aid in finding amplification factors for unusual cases, 
Fig. 3.8.2 is given to show the effect of asymmetry of an acceleration pulse. 
The pulse is triangular in shape but the time (t p ) taken to reach the peak 
value of acceleration may have any value from zero to the total duration 
(7-2) of the pulse. 

PART IV 

DISTRIBUTED MASS AND ELASTICITY 

4.1 Introduction 

It is important to be aware of the conditions under which the assumption 
of lumped parameters is permissible. In Parts I and II the cushioning 
medium was assumed to be massless, so that wave propagation (or surges) 
through it was ignored. Such surges will contribute to the acceleration 
imposed on the packaged article and we should be able to predict both the 
magnitudes and frequencies of the additional disturbances. If this is done, 
the information in Part III may be used to obtain at least a rough estimate 
of the resulting effects. In Part III itself the effects of accelerations were 
determined by studying the response of a system having only one degree of 
freedom; that is, an element of the packaged article was assumed to be a 
single mass supported by a massless spring. Every real element, of course, 
has an infinite number of degrees of freedom, so that it is important to 
discover the contribution, of the higher modes of vibration of an element, 
to the overall response. 

Both of these problems (distributed parameters of mass and elasticity in 
the cushioning medium and in an element of the packaged article) are 
studied in this part. One example of each type is considered, and the choice 
of the example in each case was influenced by considerations of expediency, 
namely that the mathematical derivations be relatively simple and lead to 
solutions for which not too lengthy computations are necessary to yield 
results that can be applied practically. At the same time, the examples 
chosen are believed to give some insight into several of the physical phe- 
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nomena involved. The treatment is by no means complete, but a more 
detailed investigation is beyond the scope of this paper. 

4.2 Effect of Distributed Mass and Elasticity of Cushioning 
on Acceleration of Packaged Article 

Referring to Fig. 4.2.1, we consider the packaged article, of mass mi , 
to be supported by distributed cushioning of mass m e and depth I. The 
cushioning may be a pad, say of rubber, in which case I is the pad thickness, 
or it may be a helical metal spring, in which case I is the coil length. The 
package is dropped vertically from a height // and, has attained a velocity v 
at the instant of contact (t = 0) of the outer container and the floor. The 
outer container is assumed to be heavy enough so that there is no rebound. 
A horizontal plane in the cushioning is located by a coordinate x measured 
from the end of the cushioning attached to the outer container. The vertical 



J (cushion) 




Fig. 4.2.1 — Packaged article of mass ma, supported on distributed cushioning of depth I 
and mass m c , depicted at the instant of first contact of the outer container (ni 3 ) and 
the floor. 

displacement of the plane .v is designated by u. The undamped motion of 
the cushioning after contact is governed by the one-dimensional wave 
equation: 



d u _ ■> d u 



(4.2.1) 



in which a is the velocity of propagation of longitudinal waves in the cushion- 
ing. If the cushioning is continuous, 



2 E 

a = — , 

P 



(4.2.2) 



where E is the modulus of elasticity and p is the density of the cushioning. 
If the cushioning is a helical spring, 



2 

a = 



m. 



(4.2.3) 



where k is the spring rate. 



(4.2.5) 
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The initial conditions of the system are 

r?i ' - -* 

L^jM 

The boundary conditions are 

[«]_o = 0, (4.2.6) 

«SL--*[SL (4 - 27) 

Equation (4.2.7) expresses the requirement that the force on the upper end 
of the cushioning must balance the inertia force of the packaged article. 
For continuous cushioning k( should be replaced by EA, where A is the 
cross-sectional area of the cushioning. 

A solution of (4.2.1) satisfying conditions (4.2.4) and (4.2.6) is 

00 

u = s. A n sin -2— sin to,,/, (4.2.8) 

where w„ is the » th root of a transcendental equation to be obtained from 

(4.2.7) and A „ is a constant to be determined by (4.2.5). Substituting 

(4.2.8) in (4.2.7) and equating coefficients of like terms of the series, we 
obtain the transcendental equation 

^tan^ = - c . (4.2.9) 

a a mi 

Substituting (4.2.8) in (4.2.5) we obtain, by the usual methods of expansion 
into trigonometric series, 

A n = - - / . ^ r~X • (4-2.10) 



w« 



(u n l . i . 2w„A" 



Hence the complete solution of the problem is 



2v sin sin u> n / 



" - ~ r / / a -n\ ■ (42 - n) 

- \— + ■ * s,n — ; 

Our chief interest is in the acceleration of w 2 . Making use of (4.2.7) 
and (4.2.9) we find, from (4.2.11), that this acceleration is 

[^l =^ i> n sina, n /, (4.2.12) 
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where 

wo = — (4.2.13) 



;//■. 



and 



2 



y m 2 V"! fl2 / 



B n = ' "" *p - » / . (4.2.14) 

inc. 1K . <a\l 

m% ml a 2 

The acceleration of m 2 is, therefore, a sum of sinusoids of frequency co„ 
and amplitude vu B n . Now, vu is the maximum acceleration that m 2 
would attain if the mass of the cushioning were negligible. Calling G n 
the maximum acceleration in the w ,h mode and G Q the maximum accelera- 
tion neglecting the mass of the cushioning, as in Part I, we have 

^ = B n . (4.2.15) 

But B n depends only on the ratio m c /m 2 , as may be seen from equations 
(4.2.9) and (4.2.14). Similarly the ratio of the frequency (co„) of any mode 
to the frequency (« ) with massless cushioning depends only on m c /m 2 , 
as may be seen from equations (4.2.3), (4.2.9) and (4.2.13). Hence, both 
the amplitude and frequency ratios for the acceleration in any mode depend 
only on the ratio of the mass of the cushioning to the mass of the packaged article. 
The ratios G n /G andco„/a>o are plotted against m c /m 2 in Figs. 4.2.2 and 4.2.3 
for the first five modes. // may be seen from these figures that the accelerations 
in the higher modes can be very important. For example, if the cushioning 
weighs half as much as the packaged article the maximum acceleration in 
the second mode is about 40% of the acceleration in the first mode and the 
latter is about the same as found by the elementary method of Part I. 
This could have a disastrous effect on an element of the packaged article if 
the latter had a fundamental frequency near that of the second mode of the 
cushioning, the latter being found, from Fig. 4.2.3, to be about five times 
the fundamental frequency of the package. 

It must be remembered that damping has been neglected in the above 
investigation and damping in the cushioning will serve to mitigate the se- 
verity of the higher mode accelerations to a great extent. However, the 
danger is always present at the start of a design and the possibilities of un- 
favorable combinations should be studied in every case. 
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Fig. 4.2.2— Influence of ratio of mass of cushioning (»»„) to mass of packaged article 
(mjj) on acceleration ratio. The numerator of the acceleration ratio is the maximum 
acceleration (G„) in the n th mode of vibration transmitted through the cu shioning. The 
denominator of the acceleration ratio is the maximum acceleration (Go = V2A*i/fKjgJ that 
the mass m 2 would experience if the mass of the cushioning were negligible. See equations 
(4.2.15), (4.2.14), (4.2.9). 




0.2 04 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



Fig. 4.2.3— Influence of ratio of mass of cushioning (m e ) to mass of packaged article 
on frequency ratio. The numerator of the frequency ratio is the frequency («„) of the 
n tb mode of vibration transmitted throug h the cushioning. The denominator of the 
freauencv ratio is the frequency (o> D = vM of vibration of the mass ms neglecting 
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4.3 Effect of Distributed Mass and Elasticity, of an Element 
of the Packaged Article, on the Amplification Factor for 
a Half-Sine -Wave Pulse Acceleration 

In this section we shall determine the contribution of the higher modes of 
vibration of a structural element to its total response to a half-sine-wave 
pulse acceleration. For the shape of the element, we choose a prismatic 
bar because this leads to the simplest mathematical formulation of the 
problem and such a bar is also a common structural element. Other con- 
siderations influence the choice of direction of acceleration with respect to 
the axis of the bar. The transverse direction (cantilever) is the most 
practical from a physical standpoint, but, for purposes of comparison with 
the one-degree-of -freedom system, the parallel (axial) direction of accelera- 
tion is the more logical. Both problems lead to solutions in the form of 
infinite series, but, in the latter case, the expression for the strain at a fixed 



I 



-v 



™2 



Fig. 4.3.1 — The system studied in Section 4.3 depicted at the instant of contact with 

the floor. 

end can be summed in terms of elementary functions without difficulty. 
Since it is necessary to determine maximum values of strain over a wide 
range of frequency ratios for the plotting of an amplification factor curve, an 
enormous reduction in the time required for accurate computations is 
obtained by choosing the axial case. Furthermore, the axial case appears to 
contain the essential features which might result in differences between the 
response of a one-degree-of-freedom system and a continuous one. 

The complete system to be studied is illustrated in Fig. 4.3.1. To the 
mass ?«2 , supported on massless cushioning of constant spring rate k 2 , 
is attached one end of an elastic prismatic bar, of length I, cross sectional 
area A , modulus of elasticity E, and density p, with its axis oriented verti- 
cally. The system is dropped from a height h so that its velocity is v 
at the instant of contact of the cushioning with the floor. The mass of the 
bar is supposed to be small in comparison with w 2 and perfect rebound is 
assumed, so that the motion of m 2 during contact is a half-sine wave of 
frequency 
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a>2 = \l — 




(4.3.1) 



The maximum acceleration of ra 2 is thus r«2 . If this magnitude of ac- 
celeration were reached very slowly, so as not to excite transient longitudinal 
waves in the bar, the maximum force between the bar and m 2 would be the 
product of the acceleration and the mass of the bar: 

F = wa 2 pAt. (4.3.2) 

Hence the strain at the end of the bar attached to m 2 would be 

VO)* pi 



€0 = 



E 



(4.3.3) 



Our problem is to find the ratio of the maximum transient strain to c . 




1 = ffa 
U 2 2 ^ 2 J 

Fig. 4.3.2— Amplification factors for an element of the packaged article having dis- 
tributed mass and elasticity. The package has linear undamped cushioning and perfect 
rebound. See equations (4.3.15), and (4.3.14). 

Let u be the displacement of a transverse plane section of the bar distant x 
from the end attached to m 2 . Then, the equation of motion of the bar is 



d u 



2 d u 

!w ~ a Hit?' 



(4.3.4) 



where a is the velocity of propagation of longitudinal waves in the bar: 

(4.3.5) 



2 E 

a = - . 



Taking the instant of first contact of the cushioning with the floor to be 
/ = 0, we know, from Part II, that the system will leave the floor when / = 
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7r/a>2 . We shall therefore treat separately, as in Part III, the motion 
during contact 

^ / 3 * 



0J 2 



and after rebound 



/ > — . 



During the first interval, the initial and boundary conditions are 

Mi_o = 0, (4.3.6) 



m 



,, , = -*> (4-3.7) 

[u] I=0 = — sin wW, (4.3.8) 

0)2 

The first and second conditions state that, at the instant of contact, all 
points in the bar are moving with the approach velocity v, without relative 
displacement. The third condition prescribes the half-sine wave motion 
of the end of the bar that is attached to m 2 . The fourth condition states 
that the strain at the free end of the bar is always zero. 

By the usual methods, a solution of (4.3.4) satisfying conditions (4.3.6) 
to (4.3.9) is found to be 

W2,« , . . . nirx . meat 

v cos — (I — x) sin o>2 / _ , „ sin — sin 



L\2C02 



^2^ 7r 2 a B-1,3.5-.- S/mra\- 1 



_ T «Jfl 



0>2 ZCdo' / 



The displacement is seen to be a forced vibration at the frequency (to 2 ) 
of the applied acceleration, on which are superposed the free vibrations of 
the bar given by the series expression. The frequency of the fundamental 
mode of vibration of the bar is wa/2 C and the frequencies of the higher modes 
are the odd integral multiples of the fundamental. 
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The strain at the attached end of the bar is 



-T-l 



= n tan 



Cd2^ . 



sin 



sin W2 



i-i E 



nirat 
IT 



™-[(5)'-']j 

/ 7T W7Tfl \ 



(4.3.11) 



W7Tfl 
2«2^ 

It may be verified that the sum of the series in (4.3.11) is given by 
nirat 



sin 



11 



'-**" {(g)'- 1 ] 



= tan — sin 0)2/ + cos co 2 t — 1, 



(4.3.12) 



(•<'<?)■ 



It should be observed that the summation is valid only in the interval < 
/ < 21/ a. However, the series is periodic with half period 21/ a and includes 
only the odd terms, so that the function repeats itself with reversed sign 
after each interval 21/ a. Hence the summation, valid for all /, can be 
written 

nirat 
4 A Shl "27 



""-•{(25) 2 - 1 ] 

= (-1)* [tan ^ sin m (t - ^\ + cos o> 2 (t - ^j- l] (4.3.13) 

! \ u 2nU ^ . . 2(m + 1)/ 
■ k = m when < 2 < 

I. \ - . a a 

> m = 0, 1, 2, 3, • • • . 

V 

We may, therefore, rewrite (4.3.11) in the form 

£ = -tan ^ sin « 2 * + (-1)* ["tan ^ sin ojt - ™) 

v a [_ a \ a / 

+ cos CO2 



('-tH 



(4.3.14) 
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0)2 

. , 2mt ^ . ^ 2(tn + 1)/ 

k — m when — < t < ! — - 

a a 

m = 0, 1, 2, 3- • - . 

The expression (4.3.14) is simple enough so that the maximum value (e ro ) 
of the strain at the attached end can be obtained without difficulty for any 
ratio of the fundamental frequency (o>i = ira/2() of the bar to the frequency 
(w 2 ) of the disturbing acceleration. The amplification factor 



A m — 



\t m \E 2ona I, ,. , ... 
j- = e m , (4.3.1a) 

VU 2 pC TT102V 



may then be calculated. The results of these calculations are plotted in 
Fig. 4.3.2. The important feature of this curve is that the amplification 
factor is everywhere less than the corresponding amplification factor for the 
one-degree-of-freedom system (Fig. 3.2.2, ft = 0). Hence the assumption 
of lumped parameters is on the side of safety as regards amplification factor. 
It is interesting to observe that the curve of A m vs. wi/« 2 , for this case, 
is a straight line between a>i/w 2 = and wi/o> 2 = 1. This arises from the 
fact that, for a>i/a> 2 ^ 1, equation (4.3.14) reduces to 



= cosw 2 / - 1, (-i ^ lj. (4.3.16) 



v 



Hence, when the duration of shock is less than the half period of the funda- 
mental mode of vibration, the maximum value of strain occurs at the end of 
impact and is equal to twice the ratio of the approach velocity to the velocity 
of wave propagation in the bar. 

The whole solution of the problem is not yet completed; for, although it is 
fairly evident from the fact that there is at least one maximum in the inter- 
val ^ t ^ t/u-2 for all values of «i/to 2 , it must be verified that the maxi- 
mum strain (and, therefore, the amplification factor) is never greater after 
/ = 7r/w 2 than before. Defining a new time coordinate 



7T 



*' = /--, (4.3.17) 



o> 2 



we have, for the initial and boundary conditions of equation (4.3.4) for 
/ > ir/«2 , 
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nair 2 , nirx 
„ sin r— sin -rr 

JY.— -Y - 1] 



'L\W/ 



v cos — (( "- .t) . » cos — ~. sin 



«C7r 2 . nirx 

(4.3.19) 



dttl a , 4 -o y^ 2co 2 l 21 



[«].-o = < (4.3.20) 

— 1 = 0. (4.3.21) 






The first and second conditions state that the displacement and velocity of 
every point in the bar must be the same at the beginning of the second inter- 
val as at the end of the first interval; the expressions in (4.3.18) and (4.3.19) 
are obtained from (4.3.10). The third condition prescribes the constant 
velocity of departure from the floor of the mass m% and, therefore, of the 
end of the bar attached to it. The last condition states, again, that the 
strain at the free end of the bar is zero. 

It may be verified that a solution of (4.3.4) satisfying conditions (4.3.18) 
to (4.3.21) is 

" nirx /nirat' \ 

u = vt' + n= E C„ sin -- sin ^— + y n ) (4 3 22) 

(f 5 0,«5 7r/a> 2 ), 



where 



nair 
Sv£ sin - - 

Cn sin 7 „ = r/ _,r -, ( 4 - 3 - 23 ) 

21 



Cn COS 7„ = 



•-•[(S)'-'] 

/ Wfl7r 2 \ 

M y 1 + cos w) 



2 r/mra V ~| 

m LVwj " l J (4-3-24) 



llira 
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Hence, the strain at the attached end of the bar is 

ntrat 
Av ^ Sin H 



["-1 

|_6:i'Jj=o tea „=i.3,5- ■• r/»7ra\ 



+iE t - 



«tt^' (4.3.25) 



7Tfl ,, = 1,3,5 ■ 



/ nira V 



The two series may be summed, as before, with the result 

— = ( — 1 ) tan — - sin o> 2 1 / — - j + cos w . 1 / — - — 1—1 

+ (-1)*' [tan^sin Ui (t> - 2 -£-') + cosw 2 (V - *£\ - l] (4.3.26) 

/ 3 7r/oj 2 , l' = t — Tr/a 1 -.. 

. 2ml ^ . ^ 2(w 4- l)/ 1 

& = m when — < * < — — , 

a a 

k = w when < t < — 

a a 

m = 0, 1, 2, 3 • • • m' = 0, 1, 2, 3 • • • . 

Once more, the expression for the strain at the attached end of the bar is 
in a form suitable for rapid calculation and it can be shown the e in equation 
(4.3.26) for t 5 ir/wo is never greater than the e in equation (4.3.14) for ^ 
I ^ t/ojo for the same &>i/«s . Hence, Fig. 4.3.2 and the conclusions follow- 
ing equations (4.3.15) and (4.3.16) need not be modified. 

Notations 

A Cross sectional area of a bar element of the packaged article. Also, a 

constant of integration. 

.lo Amplification factor when the reference acceleration is Go. Ratio of 

maximum dynamic response to the response to a slowly applied ac- 
celeration of magnitude G„g. 

A,,, Amplification factor when the reference acceleration is G m . Ratio of 

maximum dynamic response to the response to a slowly applied accelera- 
tion of magnitude G,„g. 

A n In Section 1.15, the sum of all the trapezoidal areas from xt = to Xj = 

(* 2 ),>. Also, in Section 4.2, the coefficient of the nth term of a series. 

A.4 „ The area of a trapezoid with altitude A(.v 2 )„ and sides P n -i and P n . 

a xo/l in the tension spring package. Also, in Part IV, the velocity of 

propagation of longitudinal waves. 
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B A parameter of cushioning with cubic elasticity defined in equation 

(1.5.3). Also, a constant of integration. 

B n Coefficient in the » th term of a series. 

b f/l in the tension spring package. 

C A constant of integration. 

C Coefficient of the n lh term of a series. 

c A constant defined in equation (1.7.11) 

C\ Damping coefficient of an element of a packaged article. 

c t Damping coefficient of linear cushioning. 

en The elliptic cosine function. 

do Hypothetical displacement that would result if initial spring rate were 

maintained. 

d b Maximum possible displacement of packaged article in cushioning with 

tangent elasticity. 

d m Maximum displacement of packaged article. 

d' m Value of d m when k = k' . 

d, Displacement of bi-linear cushioning at which the spring rate changes 

from ko to h. 

E Modulus of elasticity. 

e In the tension spring package the stretch of a spring when the displace- 

ment is dm- 

exp ( ) e °, where e is the Naperian base 2.718- - • 

F In section 2.7, a frictional force. 

F m In the tension spring package, the maximum force on a spring. 

f In the tension spring package, the difference between / and the distance 

between hooks of an unstretched spring. 

/, Frequency of vibration of an element of the packaged article. 

/. Frequency of vibration of the packaged article on its cushioning. 

Go Hypothetical maximum acceleration (in number of times g) that would 

result if initial spring rate were maintained. 

G. A m G m or A G , i.e. the slowly applied acceleration (in number of times g) 

that will produce the same maximum response as a transient acceleration 
of maximum value G m or Go. 

G r Maximum acceleration (in number of times g) in cushioning with fric- 

tion and spring rate kp. 

G m Absolute value of maximum acceleration of packaged article in units of 

"number of times gravitational acceleration." 

G«' Value of G m when k = k'„. 
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G„ In section 1.15, the maximum acceleration (in number of times g) ex- 

perienced by the suspended mass when dropped from a height h n . In 
Part IV, the maximum acceleration (in number of times g) of the n th 
mode of vibration. 

G T Maximum acceleration (in number of times g) after rebound. 

G. Safe value of G t . 

g Gravitational acceleration. 

h Height of drop. 

/;„ In Section 1.15, the height of fall that will cause the cushioning to dis- 

place an amount {xn) n - 

K In the tension spring package, the initial spring rate of the suspension. 

In Section 2.8, the complete elliptic integral of the first kind. 

K\, Ki, Ki The initial spring rates in the three mutually perpendicular directions 
normal to the faces of the package frame. 

k In the tension spring package, the spring rate of a spring. In Section 

2.8, the modulus of an elliptic integral. 

k, k' In Section 4.3, 0, 1, 2, 3, • • • . 

#o Initial spring rate of non-linear cushioning. 

k'o Optimum value of initial spring rate ko. 

k, Spring rate of lumped elasticity of element of packaged article. 

kz Spring rate of linear cushioning. 

ku Spring rate of bilinear cushioning after bottoming. 

kp Spring rate defined in equation (2.7.7). 

L Constant defined in equation (1.8.2). 

I In the tension spring, the projection of /, on a horizontal plane. In 

Section 4.2, length of cushioning. In Section 4.3, length of element of 
packaged article. 

/, In the tension spring package, the distance between the two support 

points of a spring when the suspended article is in the equilibrium 
position. 

M Constant denned in equation (1.8.4), equal to G„,/Go- 

m Reduced mass denned in equation (2.4.5). 

m, m' In Section 4.3, 0, 1, 2, 3, • • • . 

»i\ Lumped mass of element of packaged article. 

».'2 Lumped mass of packaged article. 

fftj Lumped mass of outer container. 

m : Mass of cushioning. 

N AP. 
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n 0,1,2,3,-... 

P Force transmitted through cushioning. 

P„ Asymptotic value of force transmissible through cushioning with hyper- 
bolic tangent elasticity. 

P„ Maximum force exerted on packaged article by cushioning. 

P n In Section 1.15, the load that produces displacement («j)». 

R Force between package and floor. 

r Coefficient of cubic term in load-displacement function for cushioning 
with cubic elasticity. 

s, t, u The direction cosines of the acceleration direction with respect to the 
normals to the faces of the package frame. 

sn The elliptic sine function. 

Ti The period of vibration of the packaged article on its cushioning. 

t Time coordinate. 

/ Time of first contact of package with floor. 

l m Time at which maximum displacement or acceleration occurs. 

t r Time at which package leaves floor on rebound. 

/, Time at which the displacement reaches d a . 

u Displacement in x direction. 

v Approach velocity. 

Wi Weight of packaged article. 

Wz Weight of outer container. 

x xi — xo; relative displacement of vi\ with respect to wj 2 . 

x -i"i - -i-2- 

x' xi - Xl. 

x' Relative displacement of Wi with respect to »w 2 at time l'. 

x In the tension spring package, the perpendicular distance from an inner 
spring support point to the nearest plane, perpendicular to the displace- 
ment direction and containing four outer spring support points. 

Xi Displacement of ntu 

xi Velocity of mi. 

Xi Acceleration of in i . 

X' Displacement of »io. 

xi Velocity of m 2 . 
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x t Acceleration of m*. 

*m»x Maximum value of x. 

{x2) n In Section 1.15, the displacement associated with the « th point. 

x.i The value x would have if the acceleration reached its maximum value 
in a very long time. 

(A;v 2 )„ In Section 1.15, equals (xo) n — (at*)»_i. 

y x 2 _x 3 . 

z x«/l (tension spring package). 

«i y> 8, f , V Phase angles. 

0i Fraction of critical damping of an element of the packaged article. 

fit Fraction of critical damping of package cushioning. 

7„ Phase angle of tt* h term of series (equation (4.3.22)). 

e Strain at attached end of element under transient conditions. 

«o Strain at attached end of element under non-transient conditions. 

t n Maximum strain at attached end of element under transient conditions. 

Angle between the displacement direction and the acceleration di- 
rection. 

x 3.14159- ••. 

p Density (mass per unit of volume) 

ro Pulse duration of a half-sine-wave acceleration. 

n Pulse duration associated with non-linear cushioning. 

m Duration of bottoming of cushioning with bi-linear elasticity. 

tp Time required to reach peak value of a triangular acceleration pulse. 

w Radian frequency defined in equation (2.4.6). 

ui Radian frequency of vibration of an element of the packaged article. 

»i' Radian frequency of vibration of damped element of packaged article. 

to- Radian frequency of vibration of the packaged article on its cushioning. 

«a' Radian frequency of vibration of the packaged article on damped 
cushioning. 

oji, A frequency defined in equation (2.10.10). 

u c A frequency defined in equation (2.8.8). 

oj„ Radian frequency of nth mode. 



